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Abstract 



In this thesis we use radio-interferometric observations to probe the intensity fluctua- 
tion power spectrum of the neutral hydrogen (H I ) 21-cm emission from the interstellar 
medium (ISM) of the external galaxies. We develop a visibility based power spectrum 
estimator to probe the power spectrum directly from the interferometric observations 
and also estimate the errors in it. A numerical simulation of synthetic observations 
is also performed to access the efficacy and limitations of this estimator. We use this 
estimator to evaluate the power spectrum of three individual spiral galaxies, a dwarf 
galaxy sample and THINGS Q spiral galaxies. In each case, the power spectrum is found 
to follow a power law Phi(^) = AW^ over a specific length scale range. The estimated 
value for the slope a ranges from ~ —1.5 to ~ —2.6 for the sample of dwarf galaxies. 
We interpret this bi-modality as arising due to 2D turbulence on length scales much 
larger than the scale-height of the galaxy disk and 3D otherwise. The power law slope 
also show a weak correlation with the surface density of star formation rate for these 
galaxies. We found for the external spiral galaxies the power spectrum is a power law up 
to a length scale of 10 to 16 kpc, indicating turbulence to be operational at these large 
length scales. We measure the scale-height of the external nearly face-on spiral galaxy 
NGC 1058 to be ~ 500 pc interpreting a change in the slope observed in it's power 
spectrum. Power spectrum of the harassed galaxy NGC 4254 is found to be affected by 
galaxy harassment. For 18 spiral galaxies derived from the THINGS sample we found 
no correlation between the power law slope and the different dynamical parameters of 
these galaxies. For most of the dwarf and spiral galaxies, the intensity fluctuation power 
spectrum is found to be a direct probe of the density fluctuation and not affected by the 
velocity fluctuations in ISM. We would attempt to understand these new observations 
in terms of physical models of the ISM in future. 



Keywords: Interstellar Medium, Turbulence, Galaxies, Scale-height, Power Spec- 
trum, visibility correlation. 



* THINGS: The HI Nearby Galaxy Survey 
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1 Introduction 



It is now well established (iBinney and Merrifield 



1998 



that all galaxies have an in- 



terstellar medium (ISM) irrespective of their morphological class. This medium carries 
10 — 15% of the baryonic mass of a galaxy and is known to contain rarefied gas (mostly 
hydrogen and helium), dust (silicate and graphite), relativistic electrons and large scale 
magnetic fields. The structure of the ISM is different in the spiral and irregular galaxies 
compared to that in ellipticals. Moreover, for most of the observed galaxies, the gaseous 
ISM extends well beyond the stellar disk. Most of this gaseous medium is hydrogen 
(70%) and helium (28%), whereas there is a trace of the heavier elements (1.5%) in- 
jected to the ISM during supernovae. Interstellar matter is distributed in discrete clouds 
of molecular and atomic gas with no preferred cloud scale. The molecular clouds have a 
temperature of about 100 — 200 K compared to the cold atomic gas with temperature 
50 — 100 K. These denser clouds occupy only 1 — 2 % of the interstellar space and con- 
stitutes half of its mass. The rest of the ISM in between these clouds consists of warm 
neutral and ionized gas. These various forms of interstellar matter are highly coupled 
and is regarded as a single dynamical entity. In the process of star formation, clumps of 
ISM driven by gravity collapses to give birth to stars, whereas, dying stars in supernovae, 
ejects a huge amount of gas, dust and metal in the ISM making it highly coupled with 
the stellar disk of the galaxy. 
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1 Introduction 



1.1 A short overview on the ISM turbulence 



Over the last few decades the structure and dynamics of the ISM has been extensively 
studied by observing it in both emission and absorption. The ISM is found to have 
rich scale invariant structures arranged in clouds of different shapes and sizes. Based on 
various ob servations and theoret i cal un derstandings it is now believed that the ISM is 



turbulent. 



Elmegreen and Scalo 



( I2OO4I ) give a revie w of the progress in observa tions of 



Scalo and Elmegreen 



ISM turbulence over the last few decades whereas 
the theoretical implications and effects of these observations. 

Pio neering work in the theory of ISM was done as early as 1951 by 



(I2OOJ) review 



von Weizsacker 



fll95ll ). He outlined a model for the ISM with a hierarchy of structures formed by the 
supersonic turbulence stirred at the large length-scales and dissipated at small scales by 
atomic viscosity. In the following year one of the f irst statistical models of a continu - 



ous and correlated gas distribution was proposed by iChandrasekhar and Miinchl (119521 ) 



who successfully applied it to the extinction fluctuations of Milky Way surface bright- 



ness. Observationally. 



absorption of starligh t by t 



structure of the ISM f Binney and Merrifield 



le Na and Ca"*" ions revealed clumpy 



19981). 



Recent observations and numerical modellings, suggest that su p ersonic turbulent flows 



play a large role in the processes of star formation (ILi et al 



20031) 



Vazquez-Semadeni et al. 



( I2OO3I ) have shown, that first the transient clumpy molecular cloud structures are gener- 
ated through turbulent fragmentation and then the individual clumps collapse to form 
proto-stars. 

Turbulence generates scale-invariant stochastic fluctuations in the velocity and den- 
sity field of the ISM. These fluctuations can be traced observing the ISM in emission 
and absorption lines of gas and metals. Observationally, statistical estimators are used 
to quantify the nature of these fluctuations as well as the turbulence that generated by 
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1.1 A short overview on the ISM turbulence 



them. Velocity fluctuations are quantified using structure function of order two whereas 
auto-correlation function or the power spectrum is used to determine the density fluc- 
tuations. Most of the knowledge of ISM turbulence is acquired using these statistical 
measures. Here we briefly discuss a few statistical tools which are widely used to measure 
and quantify the ISM turbulence. 

Auto-correlation function and power spectrum 

The stochastic density field generated by fully developed turbulence is statistically homo- 
geneous and isotropic and is believed to follow the Gaussian statistics. The statistical 
nature of such fluctuations can be completely specified by either the auto-correlation 
function or the power spectrum. The auto-correlation function ^^(| 6r \) of any homo- 
geneous and isotropic scalar field A is defined as 



Fourier transform of the auto-correlation function is known as the power spectrum, i.e, 



where k =\ k \. For both the above definitions, the angular brackets denote ensemble 
averages. For a scale-invariant fluctuation both the auto-correlation function as well as 
the power spectrum are expected to follow power laws, i.e, PA{k) ~ A;". The slope of 
these power law a is believed to be a good quantifier of the nature of the turbulence 
generating these structures. 






(1.2) 
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Structure function 

The structure function of order p for an observable A is defined as 

Spi\5r\) = {\Air)-Air + 6r)n. (1.3) 

For a scale-invariant field, Sp{Sr) are also expected to be power laws, i.e, Sp{Sr) ~ Sr^^. 
It is to be noted that the structure function of order two is directly related to the 
auto-correlation function by 5*2(1 5r |) = 2(< A'^{r) > —^{\ 5r |)). 



Delta variance 

Another common measure of a stochastic field is the delta variance defined as 

al{L) = l^j'J'^ dx [{A{r + x)-<A>) Q{x)f^ , (1.4) 
where the cutoff function is given by 



2 \ 2 I 1 if X < L/2, 

\ iiL/2<x < 3L/2. 



Power spectra measurements are typically limited by sample variance at large length- 
scales. Delta variance avoids this problem by using the cutoff function Q{x). 

It is to keep in mind that the estimators discussed above do not provide complete infor- 
mation if the stochastic field they probe has non-Gaussianity. Probes for non Gaussian 
fluctuations include structure function of higher order, higher order correlation functions, 
bispectrum, trispectrum etc. However, measuring these from the present observations is 
very difficult and is not usually addressed. 
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1.1 A short overview on the ISM turbulence 

Most of the known facts about ISM turbulence is inferred by estimating the above 
quantifiers from mostly spectroscopic observations and then comparing to existing mod- 
els of turbulence. Here we discuss how some important aspects of ISM turbulence are 
inferred from various observations. 



1.1.1 Observations an inferences 



von Hoernerl (jl95l[ ) found that the structure function of velocity fluctuations in the di- 
rection of the Orion nebulae is a power law. The slope was found to be ~ 0. 25 — 0.5 which 
indicated a Kolmogorov energy spectrum. However, using a better data 



Wilson et al. 



(119591 ) got a steeper power law with slope ~ 0.6 and inferred in favour of compressible 
turbulence. 

Apart from being scale-invariant, the density field generated by turbulence is also 
intermittent. Radio scintillation observati ons showed the clurn p y nat ure of ISM at 



very small scales down to 10^ m or lower f lLittle and Matheson 



19731). These could 



be understood using a model of magnetic turbulence generated by streaming instabil- 



ities f Wentzel 



19681 ). However, the relation of turbulence at these small length-scales 
and that observed at parsec length-scales was not very well understood. 

During 1980s, a series of obs ervation of spectral l ine- w idth at different line of sight 



Gill and HenriksenI ( 1l990l ) introduced wavelet transfor- 



revealed power law correlations, 
mation for the anal ysis of ^^CO cen troid velocities in L 1551, and the measured velocity- 



size slope was 0.7. 



HobsonI (119921 ) used a clump finding algorithm and various other 



correlation techniques for HCQ '^ and HON along M 17SW. They found scale-invariant 



fluctuations for scales < 1 pc. 



Miesch and Zweibell (119941 ) used the structure function 



to determine a slope of 0.43 ± 0.15 for the velocity-size relation. 

Comparing delta variance observations with models with and without gravity 



Ossenkopf et al 
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( I2OOII ) concluded that self gravity increases power at small scales. The observations of the 



Polaries flare by 



Ossenkopf and Mac Low! (120021 ) is compared with the models of delta 



variance for isothermal MHD turbulence. They found that the turbulence is driven from 
outside and the dissipation scale is smaller than their resolution scale. 
Observing the extinction in the Taurus region, 



Padoan et al. 



(]2002h found that Cp/Cs 



varies for p = 1 to 20 in a simi 
proposed by 



Boldvrev et al 



ar manner as in a model of supersonic turbulence 



fl2002l ). In these models the dissipation of supersonic tur- 



bulence occurs in sheets giving Cp/C: 



(IShe and Leveque 



1994h . 



3- 



Padoan et al. 



p/9 



3 assuming She-Leveque scaling 



( I2OO3I ) got a similar result observing ^'^CO emis- 



sion lines along Perseus and Taurus. 

There is a considerable amount of effort in quantifying the scale-invariant structures 
present in the ISM using fractal and multi-fractal analysis. Pe rimeter area dimension ar e 



measured for the contour map s of 10 micron dust intensity (IBazell and Desert 



CO emission f lFalgarone et al 



19881) 



I991I ) etc. The estimated value for the perimeter area 



dimension was 1.2 — 1.5 which is consistent with that of a slice from the 3D structures 



observed in the terrestrial clouds and laboratory turbulence (jSreenivasan 



19911 ). 



Chappell and Scald fl200lf ) have performed multi-fractal analysis of the column density 
maps derived from IRAS (The Infrared Astronomical Satellite) 60 fi and 100 fi images. 
They observed the region to region diversity in the multi-fractal spectra was significantly 
different from that expected from incompressible turbulence. 

These observations not only established that the ISM has clumpy scale-invariant struc- 
tures generated by turbulence over a large length-scale range, but also provided hints for 
the characteristics of the ISM turbulence. Though initially some similarity was found 
with the Kolmogorov energy spectrum, it is now believed that ISM turbulence is quite 
different. It is a compressible and supersonic fluid turbulence driven at different scales by 
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different mechanisms and induces scale-invariant multi-fractal structures in the density 
as well as the velocity field. 



1.1.2 Energy sources 



There is no wide-scale effective model for ISM turbulence which can explain all the ob- 
served features. The physical processes responsible for turbulence also differ enormously 
across the length-scales. 



At smaller scales radio scintillati on observations of intermit t ency can be we 



by sonic reflection of shock waves f llkeuchi and Spitzerl. 



mic ray streaming an d other instabi 
from larger scales etc. iLazarian et al. 



i ties (Wentze. 



1984 



19691: 



Ferriere et al 



Hall 



1 exp lained 



19881), cos- 



1980l ). energy cascades 



( 120041 ) suggested that moderate turbulence can be 



maintained in the ionized c louds by variations in the background ionizing UV radiation. 



Norman and Ferrara 



( 1l996l ) recognized that most Type II SN contribute to cluster winds 
and super-bubbles at a length-scale range of 100 — 500 pc. 

The gravitational binding energ y of a galaxy disk heats the st ellar population during 



swing-amplified shear instabilities 



triggering turbulence ( iHuber and Pfenniger 



Fuchs and von Linden 



19981 ) and also heats the gas 



2OO2I ). These gas clouds co ntiguously col- 



lapse to feed more turbulence at the smaller scales (jSemelin et al. 



19991 ). This kind of 



gravitational source for turbulence is consistent with the observed power spectra of fioc- 



culent spiral arms ( lElme green et al. 



20031 ). Since gravity acts at all the length-scales, 
the self gravity can feed turbulence at a wide range of length- scales. 

Since ISM extends well beyond the stellar disk, galactic rotation driven turbulence 
can be effect ive at the outer regions of th e galaxy's d i sk. M agneto-rotational instability. 



proposed by 



Sellwood and BalbusI ( 1l999l ): 



Kim et al 



fl2003[ ). expla ins how galactic rota- 



tion pumps energy to the gas motion and hence feeds turbulence. 



Piontek and Ostriker 
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( I2OO4I ) have considered a revised model where they assumed that the reduced dissipa- 
tion enhances energy input to turbulence by magneto-rotational instability. Rotational 
energy also generates spiral shocks where the fast-movin g inter-spiral medium hits the 



19691 ) making the ISM turbulent at 



slower dense gas in the density wave arms (iRobertsl . 
the outer edges of the disk. 

These physical mechanisms though could explain the observed fluctuations or shapes 
of the spectral lines at different length-scales, their interplay in generating wide scale 



coherent structures is not very well understood. 



1.1.3 Simulations 



Simulation play an effective role in understanding various physical models, particularly 
nonlinear processes such as turbulence. Observationally, only a static picture of ISM can 
be probed since the dynamic time scales involved are large. Using numerical hydrody- 
namic and magneto-hydrodynamic simulations it is possible to understand the effect of 
different mechanisms that drives ISM turbulence. Here we briefly discuss the different 
studies performed in this direction. 



Bania and LyonI (jl980l ) performed the first ever simulation of ISM turbulence. How- 
ever, the lack of computational facilities limited them to a very small length-scale range. 
Higher resolution studies of non-gravitating and non-rn agnetogenic compressible turbu- 



lence was first performed by 



Porter et al. 



1999 



20021 ). Their models were dominated 



by filaments, sheets and clusters in the compressible part of the flow and can explain 
the observed morphology of diffuse and small molecular clouds. 

The first nonmagnetic simulat ions which in cluded gravity were performed by[ 



Leorat et al, 



( ll990l ) on a 128^ grid, whereas 



Zhoul (jl993[ ) did the first non-gravitating compressible 



simulation to include magnetic field. These earlier works were expanded by considering 
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different aspects of ISM like magnetic and gravitational fields, realistic heating and cool- 



ing raechan isms, stellar energy injections and galactic rotation f lVazquez-Semadeni et al 



1995 



19961 ). These simulations were meant to reproduce the ISM dy namics at the 



ength - 



20031) 



scale range of ~ 1 pc. In recent times higher resolution simulations (iKaneda et al. 
observed a constant energy flux over a large length-scale range with significant departure 
from Kolmogorov energy spectrum. 

Despite the efforts simulation of ISM turbulence remained challenging due to the 
huge range of length-scales involved and high Reynolds number. It is important to 
include a wide range of length scales along with different models of energy inputs in 
the simulations to understand the interplay and cascading of the energy driving the ISM 
turbulence. Considering viscosity as the damping mechanism, Reynolds number (-Re) for 
the ISM is estimated as high as 10^ to 10^. The number of degrees of freedom for three 
dimensional incompressible turbulence scales like Re^^^. Hence simulating turbulence 
with Re ~ 10^ requires a dynamic range of 10^ which is far beyond the reach of present 
day computers. 



1.2 Radio-interferometric observation as a probe of the 
ISM turbulence 

In the previous section we have discussed the observation of ISM using different atomic 
and molecular lines which suggest turbulence. However, these elements, like dust, Ca"*", 
Na, CO etc., constitute a very small part of the interstellar matter both by mass and 
volume. These also are not very evenly distributed over the galaxy's disk. Almost 70 % 
of the gas in ISM is dominated by the atomic (H l) and ionized (H ll) hydrogen which 
is spread all over the galaxy's disk making it a better tracer of the ISM dynamics. 
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21-cm radiation 

The hyperfine structure transition radiation (21-cm radiation) of the neutral atomic 
hydrogen in the ground state (H l) is an excellent probe of the column density fluctuations 
in the ISM. This spectral line has a frequency of 1420 MHz, which falls in the radio band 
of the electromagnetic spectrum. The associated transition probability being very small, 
this line is considered as forbidden for laboratory experiments. However, at astrophysical 
scales it is possible to observe this radiation using very sensitive radio telescopes. The 
ISM has relatively small optical depth at these frequencies and the observed intensity 
fluctuations in the sky I {6, u) directly probes the column density fluctuations Nhi{0, i^) 
in ISM, i.e, 

1(6, u) = -^huA2iNHi{e)<f){u) (1.6) 

Here, is a vector in the sky, u is the observing frequency and A21 is the corresponding 
Einstein coefficient. 0(i^) gives the shape of the H I spectral line. 



Radio-interferometric observations 



Most of our understanding regarding the ISM is based on the observations in our Galaxy, 
the Milky Way. This is because the single dish radio telescopes lack the necessary 
resolution to resolve the structures in the external galaxies. Radio-interferometers like 
the VIA (Very Large Arra}!:]), GMRT (Giant Meterwave Radio Telescopcl), WSRT 
(Westerbork Synthesis Radio Telescopq^l) on the other hand provides higher resolution. 
These are capable of probing ISM structures and dynamics at smaller length-scales in 
our Galaxy and also provides means of observing turbulence at larger length-scales in 
external galaxies. In this section we briefly discuss the different efforts made in probing 



*Very Large Array, NRAO, New Mexico 
Giant Meterwave Radio Telescope, NCRA, TIFR, India 

■l-Westerbork Synthesis Radio Telescope, Netherlands Foundation for Radio Astronomy 



10 



1.2 Radio- inter ferometric observation as a probe of the ISM turbulence 



turbulence using radio-interferometric observations. 

Unlike conventional optical, infrared or ultraviolet telescopes, radio interferometers do 
not measure the sky intensity, rather they measures the visibilities which are the Fourier 
transform of the sky brightness fluctuations. These visibilities can the n be used to 



construct column density maps or the image of the observed field of view f lPerley et al 



19891 ). Using these column density maps the auto-correlation functions can be estimated. 
However, these technique has limited utility due to the the presence of correlated noise 
at the image-plane generated because of limited baseline coverages (see Appendix lAl). 



Crovisier and Dickeyl (119831 ) measured the power spectrum of H I emission from our 
galaxy squaring the directly observed visibilities (from radio-interferometric observations 
using WSRT) and subtracting the noise bias estimated from the line-free channels. They 
found that the power spectrum is roughly a power law with a slope o f ^ — 3 up to the 



GreenI ( 119931 ) estimated 



largest of 10 pc they could probe. Using similar techniques, 
the power spectrum of the H I emission of our galaxy from observations with DRAO 
(Dominion Radio Astrophysical Observator}|j) synthesis radio telescopes. He found that 
the power spectrum is a power law with slope —2.2 — 3.0 up to a length-scale of 15 pc. 



Deshpande et al. 



(120001 ) used a high dynamic range estimator to measure the power 
spectrum of the optical depth map created by observing H I in the Perseus spiral arm 
of our Galaxy in absorption towards Cassiopeia A. They found that the power spectrum 
can be well fitted by a power law with slope —2.75 for the length-scale ranges 0.10 — 3.0 
pc. Based on their observation they concluded that the ISM exhibits non-Kolmogorovian 
turbulence at these length-scales. 

Fourier transform power spectra of th e Large and Sma . 



SMC) are also found to obey power laws ( lElmegreen et al. 



1 Magellanic Clouds 



2001 



Burkhart et al. 



.MC, 



2010|) 



for a length-scale range of 20 — 200 pc. Other techniques adopted by different authors to 



^Dominion Radio Astrophysical Observatory : National Research Council, Canada 
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probe the ISM turbule nce include velocity channel an alysis or veloc i ty co ordinate spec- 



trum ( iLazarian 



correlation function f Rosolowsky et al. 



1995), fractal dimensional a nalysis (IStutzki et al. 



19981 ) and spectral 



19991 ) etc. These methods were limited either 



by the correlated noise induced by incomplete baseline coverages of radio-interferometric 
observations or by the noise bias in the power spectrum. They lack the necessary sensi- 
tivity to probe the faint emissions from the nearby dwarf or external galaxies. 



1.3 Motivations and objectives 

ISM of the Milky way is studied at length-scales ranging from 10'' m to ~ 10 pc and is 
found to have scale invariant power spectrum with slope —2.0 to —3.0 indicating tur- 
bulence is operational. For the nearby galaxies (LMC, SMC), a similar power spectrum 
is observed till a length-scale of ~ 200 pc. The proposed driving mechanisms of these 
observed turbulence in our galaxy includes galactic rotation and self gravity which oper- 
ates at even larger length-scales compared to the galaxies disk. It is necessary to probe 
the ISM of the external faint dwarf and nearby large spiral galaxies to investigate the 
existence of turbulence at large length-scales. 

Conventionally, star formation in galaxies was known to be governed by the interplay 
of the gravity and magnetic pressure modulated by the ambipolar diffusion. However, re- 



cent observations and numerical modelings suggest that supersonic 

lis process ( 



turbulent flow rather 
Star 



Li et al 



2003 



20031) 



than the static magnetic field has larger role to play in t 

formation is now believed to be a two phase process (1 Vazquez- Semadeni et al. 
first the generation of the transient clumpy molecular cloud structures by turbulent frag- 
mentation and then the collapse of the individual clumps to form the proto-stars. The 
fragmentation governed by turbulence in the ISM leads to cloud structures which are 
reflected in the fluctuations of the H I 21-cm intensity. Hence, these fragmentation can 
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be traced estimating the power spectrum of the H I 21-cm intensity fluctuations. Dwarf 
galaxies often have irregular morphologies and peculiar large scale kinematics. This also 
makes them interesting to study the ISM turbulence beside investigating the galaxy for- 
mation scenario. Observations of the nearby faint dwarf galaxies with sufficient spatial 
resolution offers us a good opportunity to understand the early stages of the star as well 
as the galaxy formation. 

Observationally it is difficult to estimate the power spectrum of the ex ternal galaxies 
becau se of the instrumental noise. A visibility based estimator was used byL 



Hobson and Maisinger 



(|2002[ ) to es timate the Cosmic Microwa ve Background Radiation (CMBR) angular power 



spectrum. iBharadwaj and Alii (120051 ) have introduced a visibility based power spec- 
trum estimator which can be used to probe the intensity fluctuations of H I 21-cm 
lin e using radio-interfero metric observations. A modified form of this estimator is used 



by 



Begum et al. 



(l2006al ) to estimate the power spectrum of the nearby dwarf galaxy 
DDO 210. The power spectrum was found to be power law with slope ~ —2.8 for a 
length-scale range till 500 pc suggesting the presence of turbulence in the ISM of the 
dwarf galaxies. 

Using a similar estimator from radio-interferometric observations of H I, it is possible 
to probe the ISM in external dwarfs and spirals and investigate the existence and nature 
of turbulence at kpc length-scales. The work presented in this thesis aims to use these 
observations to estimate the power spectrum of H I column density fluctuations in a 
large sample of galaxies. Earlier studies indicate that these fluctuations, which typically 
have a scale-invariant power spectrum, are possibly the outcome of turbulence in the 
ISM. The results of our study are expected to shed light on the nature of ISM turbulence. 

The speciflc objectives of this thesis are highlighted below: 



• Visibilities are the quantities directly measured in radio-interferometric observa- 
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tions, and there are several advantages if these be used to directly estimate the 
power spectrum. This is particularly important for the H I emission from faint 
external galaxies where the signal is weak and buried in noise. It is thus necessary 
to have a well quantified technique for visibility based power spectrum estimation. 

• Simulations play an important role in assessing the efficacy and limitations of any 
technique. It is very important to apply the power spectrum estimation technique 
to simulated galaxies in order to asses the dynamical range and errors in the 
estimate. 

• Having established a technique, it is then appropriate to apply this to estimate 
the power spectrum of the ISM of a large sample of dwarf and spiral galaxies. 

• Finally, it is desirable to have some understanding of the measured power spectra 
in terms of the underlying ISM turbulence and the processes driving it. 

1.4 Outline of the thesis 

An outline of the rest of the thesis is as follows. 

Throughout our different investigations presented in this thesis we use archival radio- 
interferometric data observed using the GMRT or the VLA. The directly observed quan- 
tity in these observations is visibility, which is the Fourier transform of the sky intensity 
fluctuations. In Chapter [2]we discuss a visibility based power spectrum estimator. The 
application of such an estimator to the H I 21-cm observation of the ISM and the error 
estimation are also discussed. 

We present the power spectrum analysis of the nearly face-on spiral galaxies NGC 628 
and NGC 1058 and a harassed galaxy NGC 4254 of the Virgo cluster in Chapter [51 In 
Chapter [4] we perform numerical simulations to study the efficacy and limitations of the 
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power spectrum estimator and access the dynamic range of length-scales for which it can 
be used. We present the result of our power spectrum analysis of a sample of 7 nearby 
external dwarf galaxies, and a sample of 18 external spiral galaxies from the THlNGq^ 
survey in Chapter O and [5] respectively. Finally, in Chapter [TJ we summarize the 
main results of the thesis and discuss about the future scopes. 



^THINGS : The HI Nearby Galaxy Survey, 



Walter et al. 



(120081) 
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2 Visibility based power spectrum 
estimator* 



2.1 Introduction 

Turbulence generates scale-invariant structures in the ISM density field. These struc- 
tures can be quantified through the power spectrum of H I 21-cm radiation. The power 
spectrum is expected to be a power law for scale-invariant fiuctuations. In this chapter we 
discuss a visibility based estimator which can be used to probe the H I intensity fiuctua- 
tion power spectrum of external faint dwarf or spiral galaxies using radio-interferometric 
observations. We also present the error estimation and limitations in this estimator. 



2.2 Visibility-visibility correlation 

The specific intensity of H I emission from a galaxy at the phase centre of the observation 
can be modeled as 

1(0, u) = w,{e) [h + si{e,u)] , (2.1) 

*A part of the work presented in this chapter is originahy pubhshed in the pa per titled "A study of 



interstchar medium of dwarf galaxies using H i power spectrum analysis" bv lDutta et al.l (|2009al ). 
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were 6 is the angle on the sky measured in radians from the center of the galaxy. We 
assume that the galaxy subtends a small angle so that 6 may be treated as a two 
dimensional (2D) planar vector on the sky. The H I specific intensity is modeled as 
the sum of a smooth component and a fluctuating component. Typically, I{6, u) is 
maximum at the center and declines with increasing \ \. We model this through a 
window function Wiy{0) which is defined so that 1^1^(0) = 1 at the center and has values 
1 > W^{0) > elsewhere. This multiplied by Ii, gives the smooth component of the 
specific intensity. For a face-on galaxy, the window function Wu{0) corresponds to the 
galaxy's radial profile. 

In our analysis we consider two different models for the window function of a galaxy 
of angular radius 6q. In the "top-hat" model it is assumed that the specific intensity has 
a constant value within a circular disk of radius 6q, and it abruptly falls to outside. 
The window function, in this Heaviside step function 



where it falls exponentially away from the center. Note that for both the models we 
have assumed the radial profile to be isotropic compared to the centre of the galaxy. 

It is also possible to use Wiy{6) to define a normalized window function Wj^{6) such 
that J dOW^{0) = 1. The second moment of 6 defined as J dO \ 9'^ \ Wi^{0) provides 
a good estimate of the angular extent of the window function. Here we have used the 
condition that this should have the same value Oq/2 for all models of the window function. 



W,{0) = 6(^0 - 0) 



(2.2) 



In the "exponential" model the window function has the form 




(2.3) 
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2.2 Visibility-visibility correlation 




Figure 2.1: W{U) for top-hat and exponential window functions for 60 = 1'. The vertical 
line marks OqK Note that the FWHM (Full Width at Half Maxima) of the 
two window functions are nearly the same and both window functions have 
a very small value for t/ > 9q^. 

to determine the width of the exponential window function in terms of ^o- Throughout 

we have assumed ^0 ^ 1 radian. 

We express the fluctuating component of the specific intensity as W^{6) 5Iy{6). Here 

SIu{6) is a stochastic fluctuation which is assumed to be statistically homogeneous 

and isotropic. The H I emission traces these fluctuations modulated by the window 

function which quantifies the large-scale H I distribution. In this thesis we use radio- 

interferometric observations to quantify the statistical properties of the fluctuations 

SI^{6). These fluctuations are believed to be the outcome of turbulence in the ISM. 



The visibility V{U, v) is the quantity that is directly observed in the radio- interferometric 
observations l It is the Fourier transform of the product of the antenna primary beam 
pattern A{6^v) and the specific intensity distribution of the galaxy. 

V{U,v) = j dO e~''^^^'^A{e,v)I{e,v). (2.4) 



^We assume the 2D approximation of the visibihty here. Interested readers are referred to the Ap- 
pendix [B] 
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Here U refers to a baseline, the antenna separation measured in units of the observing 
wavelength A. It is common practice to express U in units of kilo wavelength (kA). 
Throughout this paper we follow the usual radio-interferometric convention and express 
U in units of kilo wavelengths (kA). 

The angular extent of the galaxies that we consider here is much smaller than the 
primary beam, and it's effect may be ignored. We then have 

V{U,u) = W{U,u)I{u) + W{U,u)®6I{U,u)+Af{U,u) (2.5) 

where the tilde ~ denotes the Fourier transform of the corresponding quantity and ® 
denotes a convolution. In addition to the signal, each visibility also contains a system 
noise contribution M{U , v) which we have introduced in Eqn. (12. 5p . The noise in each 
visibility is a Gaussian random variable and the noise in the visibilities at two different 
baselines U and U is uncorrelated. We will discuss the property of the noise term 
^/{U, v) in detail in a following section. 

The Fourier transform of the normalized window functions are 



and 

W{U, v) = ^ (2.7) 

for the top-hat and exponential models respectively. Here Ji{x) is the Bessel function 
of order 1. These functions, shown in Figure 12.11 have the property that they peak 
around U = and fall off rapidly for U 9q^. This is a generic property of the window 
function, not restricted to just these two models. At baselines U ^ ^q"^ we may safely 
neglect the first term in Eqn. (12.51) whereby 



20 



2.2 Visibility-visibility correlation 



V{U, v) = W{U) ® 6I{U, v) + M{U, v) 



(2. 



We use the power spectrum of H I intensity fluctuations Phi{U) defined as 



{SI{U, u)6I*{U\ u)) = 6UU - U') Phi{U, u) 



(2.9) 



to quantify the statistical properties of the intensity fluctuations. Here, Sj){U — U') is a 
two dimensional Dirac delta function. The angular brackets denote an ensemble average 
over different realizations of the stochastic fluctuation. In reality, it is not possible to 
evaluate this ensemble average because a galaxy presents us with only a single realiza- 
tion. In practice we evaluate an angular average over different directions of U. This is 
expected to provide an estimate of the ensemble average for a statistically isotropic fluc- 
tuation. Note that, the power spectrum thus deflned can in general be a function of the 
observation frequency. However, for the observation of the ISM turbulence, the power 
spectrum is expected to remain constant within the frequency range of observation. We 
will explicitly check this assumption in a later chapter and will flnd deviation to it for 
the ISM in a particular galaxy. However, for the time being for most of the galaxies, 
it is safe to assume that the power spectrum is independent of u. We will drop the u 
dependence henceforth. 

The square of the visibilities can, in principle, be used to estimate Phi{U) 



V(t/,i/)V*(C/,z/) ) = W{U,u) (S> PHiiU) + {\AfiU,u)\^) 



(2.10) 



and t 



1993 



lis has been used in several earlier studies f Crovisier and Dickev 



1983 



Green 



Lazarian 



19951 ). This technique has limited utility to observations where the H I 



signal in each visibility exceeds the noise. This is because the noise variance in the 
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le-05 



U [kX] 



Figure 2.2: This shows the convolution (Eqn. (12.111) ) of a power law Phi{U) = U~^'^ 
with the product of two exponential window functions (Eqn. (I2.7p ) with 
^0 = 1' ((tt^o)"^ ~ IkA). The bold solid curve shows the original power law 
and the thin solid curve shows V{U, AU) for AU = 0. The dashed curves 
are for | AU \= 0.2,1.0 and 3.0 kA respectively from top to bottom. The 
direction of AU is parallel to that of U. Note that V{U, 0.2) ^ V{U, 0) and 
the correlation between two different visibilities is small for | AU \> (tt^q)^^- 

last term (| N'{U^v) p) introduces a positive bias in estimating the power spectrum. 

The noise bias can be of orders of magnitude larger than the power spectrum for the 

faint external galaxies considered here. In principle, it may be separately estimated 

using line-free channels and subtracted. In practice this is extremely difficult owing to 

uncertainties in the bandpass response which restrict the noise statistics to be estimated 

with the required accuracy. Attempts to subtract out the noise-bias have shown that 

it is not p ossible to do this at the level of accuracy required to detect the H I power 



spectrum (IBegum 



20061 ) 



The problem of noise bias can be avoided by correlating visibilities at two different 
baselines for which the noise is expected to be uncorrelated. We define the power 
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spectrum estimator as 

V{U,AU) = {V{U,u)V*{U + AU,u)) 

= J dU' W{U - U', v) W*{U + AU- U', v) Phi{U'). (2.11) 

Since W{U,i') falls off rapidly for U ^ 6^^ (Figure 12.11) . the window functions 
W{U - U', u) and W*{TJ + At/ - U\ u) in Eqn. (l2JJl have a substantial overlap 
only if |AC/| < {tiOq)^^ . Visibilities at two different baselines will be correlated only if 
|AC7| < (7r6'o)~^, and not beyond (Figure 12. 2p . In our analysis we restrict the difference 
in baselines to | AC/| < (7r6'o)"^ so that W{U + AU -U',u) ^ W{U - U') and the 
estimator V{U, AU) no longer depends on AU (Figure 12. 2p . We then use the visibility 
correlation estimator 

V{U) = {V{U,u)V*{U + AU,u)) 

= j dU' \ W{U-U',u) |2 Phi{U'). (2.12) 

The measured visibility correlation V{U) will, in general, be complex. The real part is 
the power spectrum of H I intensity fluctuations convolved with the square of the window 
function. A further simplification is possible at large baselines U ^ 6^^, provided 
I W{U, v) p decays much faster than the variations in PHi(f/). We then have 

V{U) = CPhi{U) (2.13) 

where C = J \ W{U, u) p dU is a constant. 

We use the real part of the estimator V{U) to estimate the power spectrum Phj{U). 
Our interpretation is restricted to the U range U ^ (tiOq)^^ where the convolution in 
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Figure 2.3: This shows how Um changes with a for an exponential window function with 
6^0 = 1'. Note that the discontinuities seen in the plot appear to be genuine 
features and not numerical artifacts, though the cause of these features is 
not clear at present. We find that a 4*^* order polynomial = {a + b * 

a + c* a"^ + d * + e * a^) / 9q provides a good fit to the above curve with the 
parameter values as a = 35.6, b = 104.5, c = 112.5, d = 48.9, ande = 7.7. 
We use this as a template for estimating Um in the real observations. 
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Figure 2.4: Effect of window function modifies the power spectrum differently for differ- 
ent power law exponents. 
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Eqn. (I2.12j) does not affect the shape of the power spectrum and Eqn. ( I2.13j) is a 
vahd approximation. Since for this case the estimator is independent of AU and hence 
we write it as V{U). In order to estimate the U range where this approximation is 
vahd we have numericahy evaluated Eqn. (I2.11h assuming Phi{U) to be a power law 
Phi{U) = AU"'. Figure [2.41 shows the results for both the top-hat and the exponential 
models with 9q = 1' (9^^ = 3.4 kA) and a = —1.5 and —2.5 which roughly spans the 
range of slopes usually observed for different galaxies. Using Um to denote the value (in 
kA) where the deviation from the original power law is 10%, we find that for the top-hat 
and exponential models respectively Um has values (3.1, 5.3) for a = — 1.5 and (4.8, 13.0) 
for a = —2.5. Note that Um depends on two parameters, namely a and ^o- Figure 12.31 
shows how Um changes with a for the exponential model with 6q = 1'. For other values 
of 6*0 we scale the value of Um in Figure 12.31 using Um oc 9q^. For a given power law 
index, the estimator V{U) gives a direct estimate of the power spectrum for U > Um- 



2.3 Method of analysis 

We use radio-interferometric spectroscopic multichannel data observed using GMRT or 
VLA for power spectrum analysis. In this section we discuss the method of estimating 
the power spectrum and its errors in practice. 



2.3.1 Visibility correlation 

The starting point of each analysis is the radio-interferometric spectroscopic multichan- 
nel raw uv fits data. We reduced the raw data in the usual way using standard tasks in 
classic AlPqj. For each observations, bad visibility points were edited out, after which 
the data was calibrated. The calibrated data is then inspected for H I emission and the 



■l-NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 



25 



2 Visibility Based Power Spectrum Estimator 

line-channels are identified. A continuum image is then made averaging all the line-free 
channels. In the next step, this model continuum is subtracted from the visibility data 
in the uv plane using the AIPS task UVSUB. The resulting continuum subtracted data 
was used for the subsequent analysis. Such a typical spectroscopic data has n emis- 
sion channels. To determine if the H I power spectrum changes with the width of the 
frequency channel, we combine successive channels to obtain a data set with n/N 
channels. The power spectrum is evaluated in each of this combined channels and then 
averaged over them for a particular U value. We use a range of A^ values for the power 
spectrum analysis. 

We measure V{U) by multiplying every visibility V{U, v) with all other visibilities 
V(U + AC/, v) within a disk | AC/ |< Ur, and then average over different U directions. 
We will discuss the considerations for choosing the value of Ud shortly. Note that 
correlation of a visibility with itself is excluded. To increase the signal to noise ratio we 
further average the correlations in logarithmic bins of U and over all frequency channels 
with H I emission. To choose a particular value of Ud-, we perform this analysis varying 
I AC/ I from I/tiOq to lesser values. Note that, for larger values of Ud-, we find that the 
imaginary part of the estimator becomes comparable to "P, while for smaller values 
of I AC/ I, the amplitude of decreases. For a particular observation we identify a 
value of I AC/ |, for which the ratio V jV^ > 10 and the amplitude of V is maximum. 
We use this value as Uc for the subsequent analysis. 

The measured visibility correlation estimator V{U) [Eqn. (I2.12p ] is the convolution of 
the actual H I power spectrum with a window function. We do not attempt to estimate 
this window function to deconvolve the power spectrum. Our approach is to estimate a 
baseline Um, such that for U > Um the effect of the convolution can be ignored. ViJJ) 
is proportional to the H I power spectrum [Eqn. (12.13^ ] at baselines U > Um- We 
estimate the H I power spectrum using only this range {U > Um)- 
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Assuming an exponential window function [Eqn. ( 12. 3j) ]. the value of Um depends on 
the parameters 60 and a. Note that we do not attempt to actually fit an exponential 
to the observed overall H I distribution and thereby determine 6*0. The galaxies are 
typically oriented randomly compared to the line of sight direction assuming different 
inclination angles. This offers the window functions with elliptical shapes in the sky 
plane, not circular disks as assumed in Section 12. 2[ For each of the individual galaxies 
we generate column density maps. We determine the extent of the major axis and the 
minor axis of the galaxy window function from these maps at a column density value of 
10^^ atoms cm~^. We use the smaller of the two angular extents as 6q in our analysis. 
As Um oc ^0^^, choosing the smaller value gives a conservative estimate of Um- We also 
estimate V{U) combining the visibilities of the line-free channels. The upper limit Uu 
where we have a reliable estimate of the H I power spectrum is determined by the 
requirement that the real part of V{U) estimated from the line-channels should be more 
than its imaginary part and also that estimated from the line-free channels. 



2.3.2 Error estimation 



The variance a y of the visibi. 
pendix [Ul also 



AhetaL 



mom 



ity correlation estimator V{U) is given by (see Ap- 



ex. 



V 



VjUf 



Nl 



(2.14) 



where Ng is the number of independent estimates of V{U) in a given U-hin, Nj, is 
the total number of visibility pairs in the bin and N2 is the noise variance. The first 
term in Eqn. ( ]C13[) is a contribution from the sample variance where the second term 
incorporates the noise contribution. We evaluate N2 estimating the variance in the 
visibility values. Since in a particular interferometric observation, the uv space can be 
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sparsely sampled, we grid the uv plane. We chose the grid size such that each grid 
have same area as a circle of radius Uo- We determine the fraction of grids filled with 
visibility values in each annular bin (/g)- We use ^^p^ fc as an estimate for Ng. Here, 
6U is the width of a bin at U from origin in the uv plane. With these parameters the 
variance in the visibility correlation estimator is given by 



Note that the first term under the square-root sign in the above equation dominates for 
smaller values of U, since there both U and 6U are small. The second term, however, 
dominates at large baselines where the uv plane is very sparsely sampled. 

2.3.3 Fitting a power law to the power spectrum 

We check if the estimated power spectrum follow a power law over a length-scale range 
and evaluate the slope of this power law. The procedure that we adopt to determine 
the best fit power law to the H I power spectrum is as follows. We first visually identify 
a baseline Um beyond which {U > Um) the visibility correlation estimator appears to 
be a power law in U. The largest baseline Umax till which the estimator can be used is 
limited by the requirement that the real part of V{U) estimated from the line-channels 
is greater than its imaginary part as well as the real part of the V{U) estimated from the 
line-free channels. We use the range Um ^ U < Umax to fit a power law V{U) = A f/" 
through a minimization. The best fit a obtained by this fitting procedure is used 
to get a revised estimate for Um [see Figure 12. 3j and the power law fitting is repeated 
with this new value of Um- We iterate this procedure a few times till it converges. 
To test whether the impact of the window function is actually small, we convolve the 
best fit power law with | ([/, i^) p assuming an Exponential window. The convolved 




(2.15) 
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power spectra are visually inspected to asses the deviations from the power law. The 
goodness of fit to the data was also estimated by calculating for the convolved power 
spectrum. We accept the final fit only after ensuring that the effect of the convolution 
can actually be ignored. We note the value of the best fit power law index a and the 
range of U values {Umin — Umax) for the fit. We also determine the values Rmax = D/Umin 
and Rmin = D/Umax, where D is the distance to the galaxy. These provide the range 
of length-scales at the galaxy's disk over which the power-law fit holds. We can then 
conclusively state that over those length-scales in the galaxy the power spectrum assumes 
power law and the ISM is turbulent. 

2.4 Discussions and Conclusions 

We discuss the major points of this chapter here and conclude. 

• In this chapter we have presented a method of estimating the power spectrum of 
H I intensity fluctuation from radio- interferometric observations using the visibility 
correlation. The effect of noise bias is avoided correlating the visibilities at nearby 
baselines. This method holds the key in measuring turbulence in the external 
galaxies. 

• We have discussed an operating formula for estimating the error bars in the visibil- 
ity correlation estimator we presented. This estimates of the error incorporates the 
effect of sample variance as well as the statistical errors and signifies the quality 
of the power spectrum measurement. 

• For a turbulent ISM, the power spectrum of the H I intensity fiuctuation is ex- 
pected to be a power law. We have discussed the procedure to obtain the best fit 
power law slope a from the power spectrum estimates. However, since we do not 
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exclusively evaluate the shape of the window function here, our method can not 
be used to measure the intensity of these fluctuations. 



Armed with this technique, in the following chapters we will use this estimator for 
external dwarf and spiral galaxies and investigate the relation of turbulence with some 
of their dynamical and morphological parameters. In a different investigation we have 
also used the same estimator to evaluate the power spectrum of the supernovae remnant 
and opacity fluctuation power spectrum of our galaxy. However, these wor ks are not 
releva nt in the objectiye of t he present thesis. Interested readers may look at 



(|2009|) and 



Rov et al 



Roy et al 



(120101 ) for further information. 
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Three Spiral Galaxiesl^ 



3.1 Introduction 



Power spectrum analysis is extensive 



1983 



Green 



1993 



Deshpande et al. 



y used for studying the 



2000 



Faison and Goss 



SM in our galaxy ( 



2001 



Crovisier and Dicke 



Braun and Kanekai 



20051 ). Being observations from inside the Galaxy, these are restricted to length scales 
of ~ 10 AU to 10 pc, where a power law power spectrum with a slope of ~ —2.8 
is seen. This is believed to be the outcome of ISM turbulence. The power spectra 



of LMC and SMC are also studied ( Stanimirovic et al. 



Stanimirovic and Lazarian 



1999 



Elmegreen et al. 



2001 



200ll ). where the ISM is probed at relatively larger length- 



scales (~ 200 pc). The power spectrum of the nearby dwarf galaxies are also found to 
be power laws of similar slope as observed for our Galaxy. 
Recently 



Begum et al 



(l2006al ) have used the power spectrum estimator defined in 



Chapter [2] to probe the ISM of the external dwarf galaxy DDO 210. The power 



'The works presented in this chapter are derived from the foU owing three o rigina l pubhcations: 



(1) "H I power spectrum of the spiral galaxy NGC 628" bv lDutta et al.l (|200S )^ 



(2) "The scale-height of NGC 1058 measured from its H i power spec t rum." by 

(3) "Turbulence in the harassed galaxy NGC 4254" by lOutta etall ()2010al) . 



Dutta et all (|2009b( ). 
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Galaxy 


Major 


Minor 


D 




Reff 




(') 


(') 


(M pc) 


n 




NGC 628 


15.0 


11.0 


8.0 


13 


1 


NGC 1058 


13.0 


11.0 


10.0 


11 


2 


NGC 4254 


8.0 


6.5 


16.7 


42 


3 



Table 3.1: Some relevant parameters of the 3 galaxies for which we do the power 
spectrum analysis in this cha pter. The refferences are as follows: 1 



(119931) 



Kamphuis and Briggsl (119921 ) , 2- iPetric and RupenI ( 120071 ) , 3- iPhookun et al. 



spectrum is found to be a power law indicating the presence of turbulence. Further, 
they found that the slope of the power spectrum is ~ —2.7, which is similar to what 
was observed for our Galaxy, LMC or SMC previously. Since DDO 210 is a relatively 
nearby faint dwarf galaxy, they have managed to probe a length-scale of 80 to 500 pc, 
which is also similar to what was probed for our Galaxy. 

Observing nearby spiral galaxies using radio interferometers like GMRT or VLA, it 
is possible to probe the ISM at larger length-scales. Such observations can directly 
estimate the operating length-scale range of the ISM turbulence. In this chapter we 
present the power spectrum analysis of 3 external spiral galaxies namely NGC 628, 
NGC 1058 and NGC 4254 (Table 13. ip . First two of these galaxies are nearly face on 
spirals whereas the galaxy NGC 4254 resides in a cluster environment and is believed 
to be a harassed galaxy. The objective of this chapter will be to investigate the ISM 
intensity fluctuation power spectrum at larger length-scales than it is hitherto probed. 
Moreover, the power spectrum of NGC 4254 can reveal the nature of ISM for the galaxies 
residing and interacting in a cluster environment. 
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3.2 Evidence of turbulence at large length-scales 
NGC 628 



NGC 628 (M 74 ) is a nearly face-on SA 



range 6° to 13° (IKamphuis and Briggs 



s)c s piral galaxy with an inclination angle in the 



19921 ) ■ It has a very large H I disk extending out 



Elme green et al. 



( I2OO6I ) have found a scale- 



to more than 3 times the Holmberg diameter 
free size and luminosity distribution of star forming regions in this galaxy, indicating 
turbulence to be functional here. The dist ance to this galaxy is u ncertain with previous 



estimates ranging from 6.5 Mpc to 10 Mpc. 



992 



Briggs et al. 



1980 



and 



Kamphuis and Briggs 



lave used a Hubble flow distance of ~ 10 Mpc. On the other hand, 



Sharina et al. 



19961 ) estimated a distance of 7.8 ± 0.9 Mpc from the brightest blue star in the galaxy. 



This distance estimate ra atches with an ind ependent photomet ric distance estimate by 



Sohn and Davidgd (119961 ). In a recent study 



Vinko et al. 



( I2OO4I ) inferred the distance to 



be 6.7 ± 4.5 Mpc by applying the expanding photosphere method to the hyper novae 
SN 2002ap. Here we adopt the photometric distance of 8 Mpc for NGC 628. At this 
distance 1" corresponds to 38.8 pc. 



3.2.1 A brief description of the data used 

We have used archival H I data of NGC 628 from the Very Large Array (VLA). The 
observations had been carried out on August 1''* and November 14*'^, 1993 respectively in 
the C and D conflgurations of the VLA, as a part of the AA0163 observing program. We 
performed the data analysis following the method described in Section 12.31 The multi- 
conflguration data were downloaded from the VLA archive and reduced in the usual 
... „ ... .„ ...e Mpfi Fo. ea. VLA .o„fl....o„. .ad v...... 

^NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 
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6 8 10 




01 34 45 30 15 00 33 45 30 15 

RIGHT ASCENSION (61950) 



Figure 3.1: The 38" x 36" resolution integrated H I column density map of NGC 628 
(contours) overlayed on the optical DSS image (grey scale). The contour 
levels are 0.24, 1.96, 3.68, 5.41, 7.13, 8.85, 10.57 and 12.29 xlO^o cm'^. 

points were edited out, after which the data were calibrated. The calibrated data for 
both configurations was combined using the AIPS task DBCON. The H I emission from 
NGC 628 spans 64 central channels (i.e, n = 64) of the 256 channel spectral cube (with 
channel width 1.29 km s"^). A continuum image was made using the average of all the 
line-free channels. The continuum from the galaxy was subtracted from the data in the 
uv plane using the AIPS task UVSUB. The resulting continuum subtracted data was 
used for the subsequent analysis. Figure [3. II shows a total H I column density (Moment 
0) map of NGC 628 from an image made from this data. The H I disc of the galaxy 
is nearly face-on. The angular extent of the H I distribut ion in Figure | 3.1lis roughly 
11' X 15'. Using deep VLA D array mosaic observations, 
had detected a ~ 38' x 31' faint diffuse H I envelope around NGC 628. This extended 
envelope is not detected in the current data set; the emission that we do detect is instead 
restricted to the main H I disc of NGC 628. The results we discuss below are hence also 



Kamphuis and Briggsl (jl992[ ) 
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U [KA] U [KA] 

Figure 3.2: A) Real and imaginary parts of the observed value of the H I power spectrum 
estimator V{U) for = 64 ie. all frequency channels with H I emission were 
collapsed into a single channel. The real part is also shown using 64 line-free 
channels collapsed into a single channel. B) Best fit power law to V{U). The 
channel width is varied (A^ = 1,32,64; top to bottom), 1 a error-bars are 
shown only for A^ = 64. 

relevant only to the gas in the main H I disc around NGC 628. 

3.2.2 Results and discussion 

Figure I3.2.1I A shows the real and imaginary parts of the observed value of the estimator 
V{U) for the 64 channels which have H I emission. As expected from the theoretical 
considerations mentioned earlier, the imaginary part is well suppressed compared to the 
real part. To test for a possible contribution from residual continuum, we also show the 
real part of V{U) using 64 line-free channels. This is found to be much smaller than 
the signal for the channels with H I emission. We find that a power law A U°' with 
slope a = —1.6 ± 0.2 provides a good fit to the the observed V{U) over the U range 
Umin = I.OKA to Umax = 10.0 K A (Figure 13.2. IB ) which corresponds to spatial scales 

of Rmin = 800pC to Rmax = 8 kpC. 

Both H I density fluctuations as well as spatial fluctuations in the velocity of the 
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N 


Av 


Umin 


Umax 


Rmin 


Rmax 


a 




(km s~^) 


(k A) 


(k A) 


(k pc) 


(k pc) 




1 


1.3 


1.0 


10.0 


0.8 


8.0 


-1.7±0.2 


32 


42.3 


1.0 


10.0 


0.8 


8.0 


-1.6 ±0.2 


64 


82.5 


1.0 


10.0 


0.8 


8.0 


-1.6 ±0.2 



Table 3.2: NGC 628 power spectrum are well fitted by a power law of Phi = AU°'. 

This table summarises the result for different channel widths. Note that the 
slope of the power spectrum does not changes with the width of the velocity 
channel. 

H I gas c ontribute to fluctuation s in th e H I specific intensity. Considering a turbu- 



lent ISM, iLazarian and PogosvanI (j2000l ) have shown that it is possible to disentangle 
these two contributions by studying the behavior of the H I power spectrum as the 
thickness of the frequency channel is varied. If the observed H I power spectrum is 
due to the gas velocities, the slope of the power spectrum is predicted to change as 
the frequency channel thickness is increased. We have tested this by determining the 
H I power spectrum for different values of the channel width in the range 1.29 km s~^ 
to 82.6 km s^^ (Figure 13.2. IB ) and did not find any change in the slope of the H I 
power spectrum (Table 13. 2p . As the thickest channel that we have used is consid- 
erably w ider than the typical H I v elocit y dispersion of 7 — lOkms"^ seen in spiral 



galaxies ( IShostak and van der Kruit 



19841 ). we conclude that the observed H I power 



spectru m of NGC 62 8 is pur ely due to density fluctuations. Our finding is similar to 



that of 



Begum et al. 



(j2006al ) who notic ed no change of t 



for the dwarf galaxy DDO 210. Further, 
behavior for LMC. 



Elme green et al 



le slo pe with channel width 



( 1200 ll ) also reported a similar 



Earlier studies of the Milky- Way, and a. 



DDO 210 ( 



2000 



Crovisier and Dickevl. 



Elmegreen et al. 



2001 



1983 



Green!. 



Begum et al. 



so of the dwarf galaxies LMC , SMC and 



1993 



Stanimirovic et al. 



1999 



Deshpande et al. 



2006aJ) have all found a power law H I power 
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3.2 Evidence of turbulence at large length-scales : NGC 628 



spectrum with slope ~ —3. On the contrary, we find a slope —1.6 ± 0.2 for NGC 628. 
This is a little more than one in excess of the earlier values. However, when comparing 
these values it should be noted that the earlier works have all measured the H I power 
spectrum at much smaller length-scales in the range 10 to 500 pc [MW (1 pc - 20 pc), 
SMC (4 pc - 30 pc), LMC (60 pc - 200 pc), DDO 210 (80 pc - 500 pc)] whereas the 
current measurement probes much larger lengt h-scales from 800 pc to 8.0 kpc. The typ- 



ical 



i I scale-heights within the Milky- Wav flLockman et al 



19901) and external galaxies (e.g. 



Narayan and Jog 



1984; 



Wouterloot et al 



20021 ) are well within 1.5 kpc. This 



implies that on the largest length-scales which we have probed, the turbulence is defi- 
nitely confined t o the p lane of the galaxy's disk and is therefore two dimensional (2D). 



Elmegreen et al. 



( 1200 ll ) have found that the HI power spectrum of LMC fiattens at 
large length-scales, which was interpreted as a transition from three dimensional to two 
dimensional turbulence. Based on these we conclude that the slope is different in our 
observation because it probes 2D turbulence, whereas the earlier observations were on 
length-scales smaller than the scale-height where we can expect three dimensional (3D) 
turbulence. To the best of our knowledge our results are the first observational determi- 
nation of the HI power spectrum of an external spiral galaxy at such large length-scales 
which are comparable to the radius of the galaxy's disk. 



Westpfahl et al. 



f ll999l ) have performed a fractal analysis using the perimeter-area 
dimension of intensity contours in HI images of several galaxies in the M 81 group. 
Of particular interest is the galaxy M 81, a spiral galaxy for which the perimeter-area 
dimension was found to be ~ 1.5 at a length-scale ~ 10 kpc. This observation is 
consistent with a power law power spectrum of slope —1.6 ±0.2 provided the assumption 
that the local dimension has the same value as the perimeter-area dimension is valid. 
It is difficult to probe the H I scale-height of external face-on galaxies. 



Padoan et al. 



fl200ll ) presented a method to probe the scale-height from a change in the slope of the 
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Spectral Correlation Function (SCF) 



Rosolowsky et al. 



(119991) ], and applied it to H I 



Elmegreen et al. 



(120011 ) 



data for the LMC to estimate the scale-height to be ~ 180 pc. 
suggested that one could use a change in the slope of the power spectrum of the density 
fluctuations to measure the scale-height of face on gas disks. Applying this method to H I 
data for the LMC they measured scale-height of 100 pc. To the best of our knowledge, 
prior to this work, there hava been no observational constraint on the H I scale-height 
of NGC 628. Since the scale-height is definitely less than 8 kpc and the power spectrum 
is found to have the same slope from 80 pc to 8 kpc, we conclude that the scale-height 



must be less than 800 pc. 



Kregel et al. 



(12004 ) present H I images of a large sample of 



edge on intermediate to late type spirals; from their data the ratio of the H I disk height 
to the radius of the H I disk (at a column density of 1 MqPc— 2) is ~ 0.06 ±0.015. From 
Fig. 13. II the disk of NGC 628 has a diameter of ~ 28 kpc at this column density. From 
the average thickness to radius ratio for edge on galaxies, one would expect NGC 628 
to have a scale-height of ~ 840 pc, which is consistent with our observation. 



3.3 Estimating scale-height of a face on galaxy 
NGC 1058 



NGC 1058 is an alm ost face-on, type S c 



the range 4° to 11° (jPetric and Rupen 



, spir al galaxy with an inclination angle in 



20071 ) . The di stance to this 



c ertain with previous est i mates ranging from 10.0 Mpc ( iBoroson 



galaxy is un- 



(( ISandage and Tammann 



19811 ) to 14.5 Mpc 



19741 )). Throughout this paper we adopt a distance of 10 
Mpc for NGC 1058. At this distance 1" corresponds to 48.5 pc. 

In this section we present a measurement of the power spectrum of H I intensity 
fluctuations of NGC 1058. The observed power spectrum shows a break with a steeper 
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3.3 Estimating scale-height of a face on galaxy : NGC 1058 



5 TO 15 




0241 00 4045 30 15 00 

RIGHT ASCENSION (B1950) 



Figure 3.3: The 13.6" x 13.6" resolution integrated H I column density map of NGC 1058 
(contours). The contour levels are 3., 10., 30., 50., 70., 100., 120., 130., and 
140. xlO^o cm~^ 

slope at length-scales smaller than the break. We use this to estimate the scale-height 
of NGC 1058 

3.3.1 A brief description of the data used 

We have used archival H I data of NGC 1058 from the Very Large Array (VLA). The 
observations had been carried out on 14*^ June 1993 in the C configurations of the 
VLA. The data was downloaded from the VLA archive and reduced in the usual way 
using standard tasks in classic AIP^. We adopted the data analysis procedure discussed 
in Section 12.31 for reducing the raw data and for the power spectrum analysis. The 
H I emission from NGC 1058 spans 24 central channels (i.e, channel 53 to channel 
76) of the 127 channel spectral cube. We have used only the central 16 channels (i.e, 
channel 57 to channel 72) with relatively higher H I emission. The frequency width of 

■'■NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 
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X 
CL 




u [ky 



Figure 3.4: Real (blue) and imaginary (green) parts of the observed value of the H I 
power spectrum estimator V{U) for = 16. la error-bars are shown only 
for the real part. The real part of the V{U) from line-free channels (red) is 
also shown. 



each channel corresponds to 2.58 kms~^. Figure 13.31 shows a total H I column density 
(Moment 0) map of NGC 1058 from an image made using the continuum subtracted 
data. The angular extent of the H I distribution in Figure 13.31 is measured to be 
11' X 13' at a column density of 10 ^^ atoms cm^^, which is ~ 4 times it's Holmberg 



diameter ( iPersic and Rephaeli 



20071). 



3.3.2 Results and discussion 

Figure 13.3.11 shows the real and imaginary parts of the observed value of the estimator 
V{U) for the 64 channels which have H I emission. As expected, the imaginary part and 
the real part of V{U) estimated using line-free channels is well suppressed compared to 
the real part from the line channels indicating that for these channels the observed 7^(11) 
may be directly interpreted as the H I power spectrum at U values that are considerably 
larger than (vr^o)"^ = 0.1 kA. 
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3.3 Estimating scale-height of a face on galaxy : NGC 1058 



Umin 


Umax 


Rmin 


Rmax 


N 


Av 


a 


(k A) 


(k A) 


(k pc) 


(k pc) 




(km s~^) 












1 


2.68 


-2.5 ±0.6 


6.5 


16.0 


0.6 


1.5 


8 


21.44 


-2.2 ±0.6 










16 


42.88 


-2.2 ±0.5 










1 


2.68 


-1.0 ±0.2 


1.5 


6.5 


1.5 


10.0 


8 


21.44 


-0.9 ±0.2 










16 


42.88 


-0.8 ±0.2 



Table 3.3: NGC 1058 power spectrum are well fitted by two power laws of Phi = AU°' 
in two different length scales. This table summerises the result for different 
channel width. 

The power spectrum is well fitted by two different power laws P{U) = AU", one 
with a = — 1.0 ± 0.2 for U = 1.0 kA to 6.5 kA (large length-scales) and another with 
a = — 2.5±0.6 for U = 6.5 kA to 16.0 kA (small length-scales). The results are tabulated 
in Table 13.31 We also find that the H I power spectrum of NGC 1058 does not exhibit 
a statistically significant change with increase of the channel thickness (Table 13. 3p . 
indicating that the observed power spectrum is due to H I density fluctuations. 

The slope a = — 1.0 ± 0.2 gives a good fit to the power spectrum measured at length- 



scales 1.5 — 10.0 kpc. The length-sc ale 10 kpc is definite. 



scale-heights within the Milkv- Wav flLockman et a. 



in external spiral galaxies (e.g. 



Narayan and Jog 



1984; 



y larger than the 



Wouterloot et al. 



t ypica l H I 



19901 ) and 



20021 ). Considering the H I disk of 
NGC 1058, it is reasonable to conclude that the slope a = —1.0 ± 0.2 is that of 2D 
turbulence in the plane of the galaxy's disk. On the other hand we have a slope of 
a = — 2.5 ± 0.6 at the small length scales 600 pc to 1.5 kpc. We expec t the scale- 



20041 ). The slope 



height of the galaxy's H I disk to be larger than 600 pc (iKregel et al. 
a = — 2.5 ±0.6 may thus be interpreted as being that of 3D turbulence on length-scales 



smaller than the thickness of the disk. We interpret the wavelength of 1.5 kpc where 
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we observe the transition from 3D to 2D turbulence in terms of the scale-height of the 
galaxy's H I disk. 

To our knowledge this is the first observational determination of the scale-height of a 
nearly face on spiral galaxy through its H I power spectrum. Further, the length scales 
at which this break in the power spectrum is seen is measured to an accuracy of the 
width of the bins in baseline (U) for calculating the power spectrum. Based on that 
we can also put 1 a errors to our estimate, whereby we have the scale-height 1.5 ± 0.3 



kpc. Since the disk of NGC 1058 
density of 1 Mq pc~^), following 



las a diameter of 3 5 kpc (Figure 13. 3p at the column 



Kregel et al 



( 120041 ) we expect it's scale-height to be 



1.05 ± 0.26 kpc, consistent with our result. Note that, H I scale-height is known to 
increase with the distance from the centre and our value is an average over the entire 
disk. 

It is interesting to note that a break in the power spectrum, like the one seen here for 



NGC 1058, 



( IRov et al. 



las a lso been observed in the shell-type supernovae remnant Cassiopeia A 



20091 ). The slope of the power spectrum changes from —2.22 ± 0.03 at 
1.6 — 10 kA (large length-scales) to —3.23 ±0.09 at 11 — 30kA (short length-scales). This 
change in the slope is interpreted as a transition from 2D to 3D magnetohydrodynamic 
turbulence. The transition occurs at a length-scale of 50 pc which corresponds to the 
thickness of the shell. 



3.4 Power spectrum of the harassed galaxy NGC 4254 



Moore et al. 



19961 ) is be- 



Galaxy harassment (frequent high speed galaxy encounters, 
lieved to be an important process in driving the morphological transformation of spiral 
galaxies to ellipticals inside clusters. Typically, the first encounters convert a normal spi- 



ral galaxy to a disturbed spiral with dramatic features drawn out from the dynamically 
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Figure 3.5: The 6.5" x 8.0" resolution integrated H I column density map of NGC 4254 
(contours) is overlayed with the optical image of the galaxy (density). The 
contour levels are 3., 5., 10., 20., 25., 35., 40. and 45. xlO^° atoms cm~^. 



cold gas. T^ 



le spiral galaxy NGC 4254, located in the nearby Virgo cluster, is found to 



have a tail f Minchin et al. 



2005f ) with neutral hydrogen (H I ) mass 2.2 x 10^ Mq within 
a distance of 120 kpc from the galaxy. This gaseous tail, with out any stellar count erpart, 



20071). Each 



is believed to be produced by an act of galaxy harassment (iHavnes et al. . 
act of harassment has the potential to induce a burst of star formation and to change 
the internal properties of the galaxy, including the properties of the ISM. In this section, 
we study the effect of the harassment on the large scale structure of the ISM. We use 
the power spectrum of H I intensity fluctuations to quantify these structures. 

NGC 4 254 is a lopsided spira l galaxy (morphological type SA(s)c), with an inclination 



of ~ 42° flPhookun et al. 



19931 ). The galaxy is located at a distance of 1 Mpc from the 



core of Virgo clu ster and is beli e ved t o be falling into the cluster with a relative velocity 



of 1300 kr n s ^ (Vollmer et al. 



16.7 Mpc flMei et al. 



20051 ). The distance to this galaxy is estimated to be 



20071 ): at this distance 1" corresponds to 81 pc. 
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3.4.1 A brief description of the data used 



We have used archival H I data of NGC 4254 from the Very Large Array (VLA). The 
observatio ns had been car r ied ou t on 4*^ — s*'* of March, 1992 using the C configuration of 



the VLA (jPhookun et ah 



19931 ). The data was downloaded from the VLA archive and 
reduced in the usual way using standard tasks in classic AlPq^. As usual, we adopted 
the data analysis procedure discussed in Section 12.31 for reducing the raw data and 
power spectrum analysis. The H I emission from NGC 4254 spans over 30 central 
channels i.e, (18 to 47) from 1401.7 MHz (2253 km s~^) to 1416.4 MHz (2562 km s"^) 
of the 63 channel spectral cube. A frequency width of 48.8 k Hz for each channel in the 
data cube corresponds to the velocity resolution of 10.32 km s~^. Figure 13.51 shows a 
integrated H I column density (Moment 0) map of NGC 4254 made from the continuum 
subtracted data. The angular extent of the H I distribution in Figure 13.51 is measured 
to be 6.5' X 8.0' at a column density o f 10^ ^ atoms cm~^, which is comparable to it's 



optical diameter (Ide Vaucouleurs et al. 



199l|). 



3.4.2 Results and discussion 

Figure 13.4.11 shows the real and imaginary parts of the observed value of the estimator 
V{U) evaluated using 16 channels from 27 to 42 with relatively H I emission. The 
imaginary part is well suppressed compared to the real part. We also estimate V{U) 
using the line-free channels to test for any contribution from the residual continuum. 
This is found to be much smaller than the signal (Figure 13.4.11) indicating that the 
continuum has been adequately subtracted out. 

The power law A U°', with a = — 1.7±0.2 is found to give a good fit to the H I power 
spectrum for the U range 2.0 kA to 10.0 kA (Table 13.41) . Corresponding largest length- 
^NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 
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R (kpc) 



2 




2 5 10 



Figure 3.6: Real and imaginary parts of the observed value of the H I power spectrum 
estimator for channels 27-42 with = 16. la error-bars are shown only for 
the real part. The real part of the estimator from line-free channels is also 
shown. 

scale (8.4 kpc) is definitely larger than the typical H I scale- heights of the galaxies. It 
is thus quite reasonable to conclude that the slope a = — 1.7 ± 0.2 is of 2D turbulence 
in the plane of the galaxy's disk. The fact that we do not observe the transition to 3D 
turbulence allows us to place an upper limit on the galaxy's scale-height as 2.6 kpc. 

Galaxy harassment is expected to have different effects on the inner and outer parts 
of the galaxy. While gas is stripped from the outer parts, the inner part looses angular 
momentum and gradually collapses to the center through repeated galaxy encounters. 



Data 


Channels 


Umin 


Umax 


a 






(k A) 


(k A) 




16 central 


27-42 


2.5 


10.0 


-1.7±0.2 


channels 










A 


23-30 


2.0 


10.0 


-2.0 ±0.3 


B 


31-38 


2.5 


10.0 


-1.5 ±0.2 


C 


39-46 


2.0 


10.0 


-2.0 ±0.3 



Table 3.4: Results of the power spectrum analysis of NGC 4254 for different data cubes 
spanning over different velocity ranges. 
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Figure 3.7: Integrated H I column density maps of the galaxy NGC 4254 using data 
cubes A, B and C. Note the diagonal movement of the centroid of emission 
from North east (A) to South west (C). The contour levels are 5., 8. and 12. 
xlO^° atoms cm~^. 




Figure 3.8: H I power spectrum for data cubes A, B and C plotted with arbitrary offsets 
to prevent them from overlapping. 
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Figure 3.9: Integrated H I column density maps of the galaxy NGC 628 using data cubes 
D, E and F. Note the diagonal movement of the centroid of emission from 
North east (D) to South west (F). The contour levels are 6., 12. and 18. 
xlO^° atoms cm~^. 
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Figure 3.10: Power spectra of the H I emission for the galaxy NGC 628 is shown for 
D (channels 108-119), E (channels 120-131) and F (channels 132-143) with 
arbitrary offsets to prevent them from overlapping. 
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We can selectively study different parts of NGC 4254, whose rotation axis is tilted at 42° 
to the line of sight, by considering different velocity channels. Our analysis till now has 
used only the central 16 channels, we now use the central 24 channels (23 — 46) for the 
subsequent analysis. We construct 3 different data cubes namely A, B and C containing 
channels 23 — 30, 31 — 38 and 39 — 46 respectively. We can now separately probe the 
North east. Central and South west parts of the galaxy (Figure 13. 7p using these three 
data cubes. For each data cube, we evaluate the H I power spectrum using individual 
channels and then average over the channels in the respective cubes. We estimate Um 
separately for A, B and C and the best fit power law is obtained for U > Um only. 
We find that for each data cube the H I power spectrum is well fitted by a power 
law (Figure 13.81) . the details being shown in Table 13. 4[ We find that the slope a is 
—2.0 ± 0.3 for A and C which probe the outer parts of the disk while it is —1.5 ± 0.2 
in B which probes the central region. To verify if this change in slope is generic to all 
galaxies, we also consider the H I power spectrum of the spiral galaxies NGC 628 and 
NGC 1058 which are not undergoing galaxy harassment. Figure 13.91 shows the regions 
of the spiral galaxy NGC 628 corresponding to each of the channel range 108-119 (D), 
120-131 (E) and 132-143 (F). We find that the power spectra of these three data cubes 
(Figure 13.101) all have the same slope ~ —1.6, which is also similar to the slope of 
the central part of NGC 4254. The results are similar for NGC 1058 and hence we do 
not explicitly show these here. Based on this we conclude that the difference in slope 
between the inner and outer parts of NGC 4254 is not a generic feature of the spiral 
galaxies rather is a consequence of the galaxy harassment. This suggests that not only 
does galaxy harassment affect the global morphology of the galaxy, it also affects the 
fine scale structure in the ISM as reflected by H I power spectrum. 

We note that NGC 4254 has an inclination of 42°, on the other hand NGC 628 and 
NGC 1058 are more face-on galaxies (inclination ~ 10°). In our analysis we selectively 
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study different parts of a galaxy by considering different velocity channel ranges. How- 
ever, since face-on galaxies do not have much range in radial velocity from the rotation 
curve, the spatial extent of H I in NGC 628 and NGC 1058 will not be very different for 
the three velocity channel ranges, which is unlike the case for NGC 4254. Hence, it is 
likely that we failed to find a change of slope in NGC 628 and NGC 1058 because of this 
effect of inclination. However, as seen in Figure [3.91 for NGC 628, all three data cubes 
D, E and F have a significant contribution of H I from the outer region of NGC 628, 
leading us to believe that the power spectrum has a slope of ~ 1.6 in the outer parts of 
NGC 628. This is similar to the value of slope seen in the central parts of NGC 4254 and 
significantly different from the slope in the outer parts of NGC 4254. Hence it indicates 
an impact of harassment in NGC 4254. Further analysis of spiral galaxies with large 
inclination angles would possibly be able to resolve this issue. 

We currently do not have an understanding of how galaxy harassment caused a steep- 
ening of the H I power spectrum in the outer parts of the galaxy. Theoretical modeling 
and the analysis of other Virgo cluster spiral galaxies are needed for further progress in 
this direction. 

3.5 Summary and conclusions 

We have presented the result of power spectrum analysis of the external spiral galaxies 
NGC 628, NGC 1058 and NGC 4254 in this chapter. Here we summarize the main 
points of this analysis. 

• For all the three galaxis we have studied here, the power spectrum is well fitted 
by a power law indicating the presence of turbulence in the ISM. The slope of the 
power law is found to be ~ -1.7 for NGC 628 and NGC 4254 and ~ -1.1 for 
NGC 1058 at large length scales. We have interpreted this as a result of the 2D 
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turbulence in the scale of the galaxy's disk. 

Our present observation also suggests that the spiral galaxies exhibit scale-invariant 
density fluctuations that extend to length-scales of ~ 10 kpc which is comparable 
to the diameter of the H I disk. While a large variety of possible energy sources 
like proto-stellar winds, supernovae, shocks, etc. have been pr oposed to drive 



turbulence seen at the smaller scales ( lElme green and Scalo 



2004h . it is still to be 



seen whether these are effective on length-scales as large as 10 kpc. 

We have successfully estimated the scale-height of the galaxy NGC 1058. This, 
to the best of our knowledge, is the flrst determination of the scale-height of any 
external face-on spiral galaxy. However, for the two other galaxies we could not 
detect any change in slope of the power spectrum and hence could only put limits 
to their scale-height. 

For all three galaxies in our sample we could not detect any change in the slope 
of the power spectrum with increasing thickness of the velocity channels. This 
definitely shows that for these galaxies v elocity does not noticeably rn odify the 



Lazarian and Pogosvanl (120001 ) shows a 



power spectrum. The result presented in 
possible modification of the power spectrum of the ISM intensity fiuctuation of 
our Galaxy. 

We observe that the outer part of the galaxy NGC 4254 has a different power 
law slope compared to the inner part. This difference, we propose, is because of 
galaxy harassment. We do not find any such difference in slope for the other two 
galaxies presented in this chapter. However, since this is a single observation of a 
harassed galaxy, it is not possible to conclusively report this as a generic feature 
of galaxy harassment. A systematic study of the other harassed galaxies in the 
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nearby Coma or Virgo cluster may provide us more interesting information about 
the ISM turbulence in such environments. 

Since it is now evident that the external spiral galaxies also have turbulent ISM, we 
may now proceed in estimating the power spectrum of a spiral galaxy sample. This 
will help us understand the nature of turbulence, its origin and the energy input to it. 
In Chapter [6] we will consider such a sample and investigate how the turbulence is 
correlated with other dynamical and morphological properties of the spiral galaxies. 



51 



3 Power Spectrum of Three Spiral Galaxies 



52 



4 Simulating Galaxy Observationsl^ 



We have presented a visibility based power spectrum estimator in Chapter [2] which we 
have used to perform power spectrum analysis of three external spiral galaxies and the 
results are presented in Chapter [31 The observed bi-modality in the value of the best 
fit power law slope a of these spiral galaxies with the same from Milkyway or nearby 
galaxies have been explained in terms of 3D and 2D turbulence at small and large scales 
respectively. The power spectrum of the galaxy NGC 1058 is found to be a broken power 
law showing a transition from 3D to 2D turbulence. The break was seen at a wavelength 
of 1.5 kpc which we interpreted as the scale-height of the galaxy. However, the actual 
length-scales associated with the scale-height could not be definitively estimated. 

In this chapter we perform numerical simulations to asses the impact of the overall 
galaxy structure on our estimates of the power spectrum of H I intensity fluctuations. 
We are particularly interested in the effect of galaxy's radial and transverse profile, 
inclination angle and the large scale density field in the power spectrum estimator we 
have been using. 



*The work presented in this chapter is originally published in the pa per titled "A study of interstellar 



medium of dwarf galaxies using H i power spectrum analysis" by IDutta et al.l (|2009aD 
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4 Simulating Galaxy Observations 

4.1 Simulation 

The starting point of our simulation is a three dimensional (3D) 512^ mesh, labeled using 
coordinates r, on which we generate a statistically homogeneous and isotropic Gaussian 
random field h{r) with power spectrum P{k) = Ck'^. We have used 7 = —2.5 throughout 
our simulations, the results can be easily generalized to other 7 values. The 3D power 
spectrum is shown in Figure 14. 1[ For reference we also show the two dimensional 
(2D) power spectrum of h{r) evaluated on a 512^ planar section of the cubic mesh. As 
expected, the 2D power spectrum also is a power law with P{k) oc k'^'^^. In all cases we 
have generated five independent realizations of the Gaussian random field, and averaged 
the power spectrum over these to reduce the statistical uncertainties. 

In our simulations we use the Gaussian random field h{r) as a model for the 3D H I 
density fiuctuations that would arise from homogeneous and isotropic 3D turbulence. 
The h{r) values on a 2D section through the cube serves as a model for the H I density 
fiuctuations in the limiting situation where we ignore the thickness of the galaxy and 
treat it as a 2D disk. Note that the resulting h{r) is a statistically homogeneous and 
isotropic Gaussian random field on the 2D section. We use this to represent the density 
fiuctuations that would arise from homogeneous and isotropic 2D turbulence. 

We next embed a galaxy in the middle of the 3D cube. The overall, large scale 3D 
structure of the galaxy is introduced through a function G{r) so that the H I density at 
any position is G{r)[hQ + h{r)]. Here G{r)ho is the smoothly varying component of the 
galaxy's H I density and G{r)h{r) is its fiuctuating component. It is assumed that the 
observer's line of sight is along the z axis. The H I density G{r)[hQ + /i(r)] is projected 
on the x — y plane. We interpret the projected values as the H I specific intensity J(0, v) 
in the plane of the sky [Eqn. ( 12.11) ]. The 512^ mesh in the x — y plane of our simulation 
corresponds to 4' x 4' on the sky. 
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4.1 Simulation 



To understand the effect of different parameters such as the shape of the window 
function, the inchnation angle etc. on the power spectrum estimation we consider two 
models for the galaxy's H I density. 



4.1.1 Case I : Without a radial profile 



We first consider a face-on galaxy, and use G{r) = exp(— 2^/2;^). This incorporates only 
the finite thickness of the disk which is characterized by the scale-height and ignores 
the galaxy's radial profile in the plane of the disk. Note that, for z < Zh this function 
closely matches the function sechj—z '^/z'^) which is also used to model the scale-height 



profile ( iBinney and Merrifield 



19981 ). For the scale-height, we have used the values 



Zh = 4: and 32 mesh units which we refer to as the "thin' and "thick" disk respectively. 



Results 

Figure 14.11 shows the power spectrum of the simulated H I specific intensity I{6, u) 
for both the thin and thick disk. We find that the H I power spectrum of the thin 
disk has a slope —2.5, same as the 3D power spectrum, at large U. There is a break 
at f/ ~ 35 kA which corresponds to l/nzh, and the slope is ~ —1.9 at smaller U. We 
interpret this change in the slope in terms of a transition from 3D fluctuations on scales 
smaller than the disk thickness (large U) to 2D fluctuations on scales larger than disk 
thickness (small U). We expect the slope to approach —1.5, the 2D slope, at very small 
U. Note that in our simulation it is not possible to evaluate the slope at very small U 
values where the sample variance is rather large. We find that the H I power spectrum 
of the thick disk is well fit by a single power law Phi{U) = AU""^'^ which has the 3D 
slope. In this case we expect the break corresponding to the transition from 3D to 2D 
at f/ ~ 4kA. This break lies in the sample variance dominated region which explains 
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U [k A] 

le+07 
le+05 

P(k) 

le+03 
10 

1 10 100 

k 

Figure 4.1: The 3D and 2D power spectrum P{k). The simulated H I power spectrum 
Phi\{U) for the thin (A) and thick (B) disk without the radial profile are also 
shown. The P{k) and k values (left and bottom axes) have been arbitrarily 
scaled to match the Phi{U) and U axes (top and right). The power-law 
fits are shown by solid lines. The different curves have been plotted with 
arbitrary offsets to make them distinguishable. 

why we do not detect it. 

4.1.2 Case II : With a radial profile 

We next incorporate the galaxy's radial profile using 

expi-z'/zl), (4.1) 

where 6 = a/x^ + is the radial coordinate in the plane of the disk. Here 9 coincides 
with the angle in the sky measured from the center of the galaxy. We have used Oq = 1' 
which corresponds to 128 mesh units . Note that the radial profile used in the simulation 
is exactly the same as the window function W{0) of the exponential model introduced in 
Chapter [21 The relative amplitude of ho and h{r) is a free parameter which decides the 
respective contributions from the smooth and the fiuctuating H I components. Using 




G{r) = exp 
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Figure 4.2: Tlie simulated H I power spectrum for tlie face-on tliick disk witli (B) no 
radial profile and (C), (D) with radial profile using 6 = 0.35, 0.7 respectively. 
(E) is same as (C) with the disk tilted at 60°. The different curves have been 
plotted with arbitrary offsets to make them distinguishable. 

the ratio R = \J {h'^)/hQ of the root mean square value of h{r) to ho to quantify this, we 
have carried out simulations for R = 0.35 and 0.7. We have also carried out simulations 
where the disk is tilted, and the normal to the disk make an angle of 60° to the line of 
sight. 

Results 

The simulated H I power spectra are shown in Figures 14.21 and 14.31 for the thick and 
thin disks respectively. We find that for a face-on disk with Oq = 1', for both the thick 
and thin disks, the results of our simulations are in good agreement with our analytic 
estimate presented in the Section 12. 2[ which predicts that the effect of the galaxy's 
radial profile is contained within a limited baseline range U < Um, which comes out to be 
U < Um. ~ 3.5 9q^, for a = —2.5 (Figure 12. 3p . The H I power spectrum is insensitive 
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Figure 4.3: The simulated H I power spectrum for the face-on thin disk with (A) no 
radial profile and (F), (G) with radial profile using 6 = 0.35, 0.7 respectively. 
(H) is same as (F) with the disk tilted at 60°. The different curves have been 
plotted with arbitrary offsets to make them distinguishable. 

to the galaxy's radial profile at f/ > Um, and the shape of the power spectrum is the 
same whether we include the galaxy's radial profile or not. We find that the value of 
Um is not very sensitive to changes in R, the ratio of the fluctuating component to the 
smooth component. For a disk of the same size, the value of Um increases when the disk 
is tilted. The image of the disk tilted by 60° has an anisotropic window function W{0) 
with angular radius 1' and 0.5' along the major and minor axes respectively. We find 
that for a tilted disk, the value of Um. is determined by the smaller of the two angular 
diameters, 0.5' in this case. 



4.2 Scale-height estimation 

In Chapter [3] we have observed a break in the power spectrum of the galaxy NGC 1058 
at a baseline of 6.5 k A (or at a wavelength of 1.5 kpc) and interpreted it as arising 
due to a transition from 3D to 2D turbulence. The results of our simulation suggests 
that a break in the power spectrum at a baseline U corresponds to an angular scale of 
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Zh = 1/itU and a length-scale of D/ttU at the plane of the sky {D is the distance to the 
galaxy). Following this for an adopted distance of 10 Mpc to the galaxy NGC 1058, its 
scale-height is 490 ± 90 pc. We have also observed that the galaxy NGC 628 exhibit 2D 
turbulence up to a baseline of 1.0 k A with a single power law of slope ~ —1.6. This 
imposes an upper limit of 320 pc on its scale-he ight. Note that the scale-height of 320 pc 



Dutta et al. 



(120081) because of the extra 



for NGC 628 differs from what is presented in 
factor of vr indicated by our simulations. We do not attempt to set a bound for the 
scale-height of the galaxy NGC 4254 here, since it has a considerably large inclination 
angle. Finally, we note that we have not considered the fact that the scale-height varies 
with the distance from the centre of the galaxy's disk. Hence our estimate gives an 
average value of the disk scale- height. 



4.3 Summary and conclusions 

We discuss the main findings of this chapter here. 

• We find that if the H I distribution is a thin disk, the power spectrum is well fitted 
using a broken power law with the change in scope at a length-scale corresponding 
to the line of sight thickness of the disk. We interpret this as a geometrical effect 
arising because of a transition from 3D to 2D fluctuation. Since the spiral galaxies 
have disk like geometry, a break in the power spectrum at a t/ value greater than 
Um is indicative of a transition from 3D to 2D at the length-scale corresponding 
to the scale-height Zh = D/nlf. 

• Our simulations demonstrate that for baselines U > Um the galaxy's radial profile 
has no effect on the power spectrum of H I fluctuations. In other words, the power 
spectrum would be unchanged if the same fluctuations were present in an uniform 
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disk with no radial profile. 

• We find that the inclination angle of the Galaxy further restricts the usable range 
of baselines for the power spectrum estimation. Our simulation shows that the 
semi-minor axis of the projected image of the galaxy has to be used to determine 
the Um value over which the power spectrum can be estimated. 

Finally, we note that we have ignored velocity fluctuations and the galaxy's rotation 
throughout our simulations. This would rearrange the H I emission amongst the different 
frequency channels. Our simulation corresponds to the situation where there is just 
a single frequency channel whose width encompasses the entire HI emission, and the 
velocities have no effect in this situation. 
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5 Power Spectrum as a Probe of 



Dwarf Galaxy Propertied^ 



5.1 Introduction 

This chapter presents the results of our analysis of a sample of 7 nearby gas rich faint 
{Mb > —13.0 mag) dwarfs. The corresponding H I radio-interferometric data used here 
is from GMRT. The observation parameters and other details are given in the Table [^^ 



Note that the galaxy GR 8 is known to have peculiar kinematics ( iBegum and Chengalui 



2003 ) whereas the spiral dwarf galaxy NGC 3741 has an extended H I disk f lBegum et al 



20051 ). These were included in our sample to see the effect of the peculiar dynamics and 



morphology of the ISM in the power spectrum. 



5.2 Notes on the individual galaxies 



We briefly discuss some of the relevant properties of the galaxies in our sample. A few 
of the galaxy parameters are also summarized in Table 15.21 which contains - Column 



*The work presented in this chapter is originally published in the pa per titled "A study of interstellar 
medium of dwarf galaxies using H i power spectrum analysis" by iDutta et al. ( 2009a ) . 
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Galaxy 


a(J2000) 


(5(J2000) 


Date of 


Time on 


References 




(h m s) 


(°"') 


observation 


source 




DDO 210 


20 46 53.00 


-12 50 57 


13-15 July 02 


16 hrs 


2,5 


NGC 3741 


1136 06.40 


+451707 


22 Jul & 26 Aug 04 


8 hrs 


3,5 


UGC 4459 


08 34 06.50 


+6610 45 


15,23,24 Nov 02 


14 hrs 


4,5 


GR 8 


12 58 40.40 


+1413 03 


16-18 Nov 02 


16 hrs 


1,5 


AND IV 


00 42 32.30 


+40 34 19 


1 Jan, 07 




6 


KK 230 


14 0710.70 


+35 03 37 


May 8, 6 Jun, 26 Nov 01 


18 hrs 


4,5 


KDG 52 


08 32 56.00 


+710146 


21-23, 27 Jun 02 


18 hrs 


4,5 



Table 5.1: Some obser vation parameters of th e galax ies i n our sample. Reference s are a s 
fol lows: 1- iBegum and Chengalurl (120031) . 2- iBegum and Chengalurl fl2004[) . 
3- iBegum et all fl2005h . 4- lBegum et al.rfl2006bh . 5- lBegum et all fl2008ah . 6- 
Chengalur et al. 2009 (in preparation) 



(1): Name of the galaxy; (2) and (3): major and minor axis in arc min calculated at a 
column density of 10^^ atoms cm~^; (4): distance in Mpc; (5): H I inclination angle; 
(6): Log[SFR] in MQyr"^ measured from Ha emission; (7): Hi velocity dispersion in 
km s~^; (8): total dynamical mass in 10^ Mq] (9): references. 

DDO 210 

DDO 210 is th e faintest [Mb ~ —10.9) relatively close (at a distance 950 ± 50 kpc; 



Lee et al 



19991 ) gas-rich member of the Local Group. The H I disc of the galaxy is 
nearly face-on. On large scales, the H I distribution is not axisymmetric; the integrated 
H I column density contours are elongated towards the east and south. No Ha emission 
was detected indicating a lack of on-going star formation in this galaxy. 



NGC 3741 



NGC 3741 is a nearby dwarf irregular gala xy (Mr ^ —13.13 ) with a gas disk that 



extends to ~ 8.8 times the Holmberg radius fiBegum et al. 



2005 



2008b|). The galaxy is 
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Galaxy 


Major 


Minor 


D 




log(SFR) 






Ref. 




(') 


(') 


(M pc) 


n 


(Mo yr-i) 


(km s~^) 


(108 Mq) 




DDO 210 


5.0 


4.5 


1.0 


26.0 


> -5.4 


6.5(1.0) 


0.34 


2,5 


NGC 3741 


27.0 


12.0 


3.0 


68.0 


-2.47 


~ 8.0 


40.3 


3,5 


UGC 4459 


4.5 


3.6 


2.1 


30.0 


-2.04 


9.0(1.6) 


1.40 


4,5 


GR 8 


4.2 


4.0 


2.1 


27.0 


-2.46 


9.0(0.8) 


0.78 


1,5 


AND IV 


15.0 


8.5 


6.7 


55.0 


-3.0 


~ 8.0 


32.5 


5,6 


KK 230 


3.0 


2.0 


1.9 


50.0 


> -5.53 


7.5(0.5) 




4,5 


KDG 52 


3.5 


3.0 


3.55 


23.0 


> -5.1 


9.0(1.0) 




4,5 



Table 5.2: Some relevant morphological and dynamical parameters of the galaxies i n 
ou r sample. The reference s are a s fol lows: 1- Begurn and Chenga lurl (2003}, 



2- 



lengg 

Begum and Chengalurl fl2004f) . 3- iBegum et al.l fl2005l ). 4- iBcgum et al 



(l2006bf ). 5- iBegum et al.l ( l2008al ). 6- Chengalur et al. 2009 (in preparation) 
fairly edge-on with kinematics inclination varying from ~ 58° — 70°. NGC 3741 appears 
to have a H I bar and is dark matter dominated with a dark to luminous mass ratio 
of ~ 149. Further, this galaxy is undergoing significant star formation in the center. 



An interp 
detail by 



ay between the neut ral ISM and star formation in this galaxy is studied in 



Begum et al. 



~ 50 kms ^ ( Gentile et al. 
motion of ~ 5 — 13kms~^ 



fl2008br). The rotation curve fiattens beyond 300" to a value 



20071). They find that the ISM in this galaxy shows radial 



Karachentsev et al. 



(l2004f ) have estimated the distance to 



this galaxy as 3.0 ± 0.3 Mpc using the tip of the red giant branch (TRGB) method. 



UGC 4459 



The faint dwarf {Mb ~ —13.37) galaxy UGC 4459 is a member of the M 81 group of 
galaxies. It is fairly isolated from its nearest neighbor UGC 4483 at a projected distance 
of 3.6° (~ 223 kpc) and a velocity difference 1 35 km s~^. UGC 4459 is a relatively 



metal poor galaxy, with 12 + log(0/H) ~ 7.62 f lKunth and Ostlin 



20001 ). The optical 



appearance of this galaxy is dominated by bright blue clumps, which emits copious 
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amount of Ha, indicating high star formation. The velocity field of UGC 4459 shows a 
large scale gradient across the galaxy with an average of ~ 4.5 km s~^ kpc~^, though 
this gradient is not consistent with that ex pected from a systematic r otating disk. This 



galaxy has a TRGB distance of 3.56 Mpc ( iKarachentsev et al. 



20041 ). 



GR 8 



'his is a faint (Mr 



( iBegum and Chengalur 



12.1) dwarf irregular galaxy with very unusual H I kinematics 



20031 ). The H I distribution in the galaxy is very clumpy, and 



shows substantial diffuse, extended gas. The high density H I clumps in the galaxy are 
associated with optical knots. In optical, it has a patchy appearance wit h the emissio n 



dominated by bright blue knots which are sites of active star formation ( iHodge 



1963). 



Both radial and circular motions are present in this galaxy. It also possesses a faint 
extended emission in Ha. The distance to this galaxy is estimated to be 2.10 Mpc 



( IKarachentsev et al. 



20041 ). 



AND IV 



This is a dwarf irregular gala xy with a moderate s urface brightness (/iv ~ 24) and a 



very blue color {V — I < 0.6 



from the center of M 31 (IMerrett et al. 



20001). It is at a projected distance of 40' 



20061 ) and has a very low ongoing star formation 



rate (SFR) of ~ O.OOlMQyr ^. Its H I disk extends to ~ 6 times the Holmberg diameter 
and also shows large scale, purely gaseous spiral arms [Chengalur et al. 2008 (in prepa- 



ratio n)]. The distance to this galaxy is estimated to be 6.7 ± 1.5 Mpc ( iFerguson et al. 



20001 ). 
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KK 230 



KK 230, the faintest(MB ~ —9.55) dwarf irregular galaxy in our samp 



periphery of the Canes Venatici I cloud of galaxies f lKarachentsev et al. 



e , lies at the 



20041 ) ■ The 



velocity field shows a gradient in the direction roughly p erpendicular to the H I and 



2006M). There 



optical major axis with a magnitude of ~ 6 km s ^ kpc ^ fiBegum et al. . 
is no measurable ongoing st ar formation in this galaxy as inferred from the absence of 



any detectable Ha emission. 



Karachentsev et al. 



( I2OO4J ) has found a tidal index of —1.0 



indicating the galaxy to be fairly isolated. They have estimated the TRGB distance to 
be 1.9 Mpc. 



KDG 52 

KDG 52 (also called M 81DwA), another member of the M81 group, is a faint dwarf 
galaxy {Mb ~ — 11.49) with a clumpy H I distribution in a broken ring surrounding the 



optical emission (IBegum et al. 



2006bl ). The H I hole is not exactly centered around the 



optical emission. This galaxy does not have any dete ctable ongoing star 



distance to this galaxy is estimated to be 3.55 Mpc (IKarachentsev et al. 
from the TRGB method. 



brmation. The 



2004h derived 



5.3 Results and discussions 

We adopt the data analysis procedure discussed in Section l2.3[ All the data are reduced 
in the usual way using standard tasks in classic AlPSB For each galaxy, after calibration, 
the frequency channels with H I emission are identified and a continuum image is made 
by combining the line-free channels. The continuum is hence subtracted from the data 



^NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 
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in the uv plane using the AIPS task UVSUB. The number of channels with H I emission 
(n) is different for each galaxy. To determine if the H I power spectrum changes with 
the width of the frequency channel, successive channels are combined to obtain n/N 
channels. Power spectrum analysis is performed for different values of in the range 
1 < iV < n . The analysis is initially carried out for = 1 and unless mentioned 
otherwise the results refer to this value. 

Figure [5. II and [ 5.21 show the results of our analysis. The results are also summarized 
in Table 16.21 We have detected the H I power spectrum of the galaxies DDO 210, 
NGC 3741, UGC 4459, GR 8 and AND IV. For these galaxies we find a range of baselines 
U, where the real part of V{U) estimated from the channels with H I emission is larger 
than the imaginary part estimated from the same channels as well as the real part 
estimated from the line-free channels. This is not true for KK 230 and KDG 52 (and 
also for AND IV with N = n), for these galaxies all three curves (Figure 15. 2p lie 
within the 1 a errors bars with very small offset and the interpretation is not straight- 
forward. For the subsequent analysis we use only the real part of V{U) estimated from 
the channels with H I emission and interpret it as power spectrum. 

We find that a power law V{U) = A U°' provides a good fit at a baseline range Umin to 
Umax with a reasonable xVi^ for DDO 210, NGC 3741, UGC 4459, GR 8 and AND IV. 
The presence of a scale-invariant, power law power spectrum indicates that turbulence 
is operational in the ISM of these galaxies. The range of length-scales for the power law 
fit differs from galaxy to galaxy and in total it covers 100 pc to 6.2 kpc. The details are 
summarized in Table 16.21 

Both H I density fluctuations as well as spatial fluctuations in the velocity of the 
H I gas contribute to fluctuations in the H I specific intensity. To test if the power 
spectrum has any contribution from the gas velocity fluctuations, we have repeated the 
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R (kpc) 




U [kX] 



Figure 5.1: Power spectrum of the galaxies DDO 210, NGC 3741 and UGC 4459. The 
real and imaginary parts of V{U) estimated after averaging N channels with 
H I emission, and the real part from N line-free channels are shown together 
for = 1 (left panel) and N = n (right panel). The error-bars are for the 
real part from channels with H I emission. The best fit power law is shown 
in bold. In each case Um is marked with a bold-faced arrow and the fit is 
restricted to U > Um where the effect of the convolution with the window 
function can be ignored. 
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Figure 5.2: Power spectrum of the galaxies GR 8, AND IV, KK 230 and M81DWA. The 
real and imaginary parts of V{U) estimated after averaging channels with 
H I emission, and the real part from A^ line-free channels are shown together 
for A^ = 1 (left panel) and N = n (right panel). The error-bars are for the 
real part from channels with H I emission. The best fit power law is shown 
in bold only where such a fit is possible. For those, Um is marked with a 
bold-faced arrow and the fit is restricted to U > Um where the effect of the 
convolution with the window function can be ignored. 
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Galaxy 


N 


Av 


Umin 


Umax 


Rmin 


Rmax 


a 


C2 






(km s"^) 


(k A) 


(k A) 


(k pc) 


(k pc) 






DDO 210 


1 


1.65 


3.7 


13.0 


0.06 


0.27 


-2.3 ±0.6 


0.0 ± 1.2 


NGC 3741 


1 


1.65 


1.6 


16.0 


0.19 


1.90 


-2.2 ±0.4 


0.0 ±0.8 


UGC 4459 


1 


1.65 


3.1 


22.0 


0.18 


1.16 


-1.8 ±0.6 


0.0 ± 1.2 


GR 8 


1 


1.65 


1.6 


23.0 


0.10 


1.30 


-1.1 ±0.4 


0.0 ±0.8 


AND IV 


1 


1.65 


0.6 


6.7 


0.56 


6.20 


-1.3±0.3 


0.0 ±0.6 


DDO 210 


10 


16.5 


3.7 


13.0 


0.06 


0.27 


-2.1 ±0.6 


0.0 ± 1.2 


NGC 3741 


8 


13.2 


1.6 


16.0 


0.19 


1.90 


-2.5 ±0.4 


0.0 ±0.8 


UGC 4459 


12 


19.8 


3.1 


22.0 


0.18 


1.16 


-1.7 ±0.4 


0.0 ± 1.2 


GR 8 


8 


13.2 


1.6 


23.0 


0.10 


1.30 


-0.7 ±0.3 


0.0 ±0.8 


AND IV 


16 


26.4 












0.0 ±0.6 



Table 5.3: The results for the 5 galaxies in our sample. Columns 1-9 give (1) Name 
of the galaxy, (2) A^: the number of channels averaged over, (3) A^: the 
corresponding velocity width, (4) and (5) Umin and Umax- the baseline range 
for which the power spectrum is determined, (6) and (7) Rmin and Rmax- 
corresponding range of length scales in the galaxies disk, (8) a: the best fit 
slope for the HI power spectrum, (7) (2- possible limits for the spectral slope 
of the velocity structure function. 
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power spectrum analysis increasing the channel thickness from = 1 to = ra. 
In addition to = 1, Figures 15.11 and 15. 2[ and Table 16.21 also show the results for 
N = n where the channel thickness spans the entire frequency range that has significant 
H I emission. We do not find a significant change in the slope of the power spectrum for 
any of the galaxies. Since for all the galaxies the thickest channel is considerably larger 
than the velocity dispersion, we conclude that the H I power spectrum is purely due 
to density fluctuations and not gas velocities. The fact that the slope does not change 
with channel thickness c an be used to constrain the v alue of the slope of the velocity 



structure function ^2 (see iLazarian and Pogosyanll2000l ). The C2 values are tabulated in 
Table lOl 

The galaxies in our sample have slope a ranging from —2.6 to —1.1. The two galaxies 

DDO 210 and NGC 3741 have slope ~ -2.5, while the slope is 1.5 for UGC 4459 

and AND IV, and —1.1 for GR 8. We have proposed a possible explanation for this 
dichotomy in the values of a in Section 12.31 This was based on the fact that DDO 210, 
where the power spectrum was measured across length-scales 100 — 500 pc, had a slope of 
—2.6 while NGC 628, a nearly face-on galaxy where the power spectrum was measured 
across length-scales 0.8 — 8 kpc had a slope of —1.6. We have interpreted the former 
as three dimensional (3D) turbulence operational at small scales whereas the latter 
was interpreted as two dimensional (2D) turbulence in the plane of the galactic disk. 
In Chapter [4] we have performed numerical simulations which tend to support this 
interpretation. Continuing with this implies that we have also measured 3D turbulence 
in NGC 3741, and 2D turbulence in UGC 4459, GR 8 and AND IV. Our dwarf galaxy 
sample contains 3 nearly face-on galaxies with ini < 30°, namely DDO 210, UGC 4459 
and GR 8. DDO 210 has 3D turbulence across the baselines 2.0 — 10.0 kA. The absence 
of a break in the power spectrum places a lower limit of 160 pc on the scale-height. The 
galaxies UGC 4459 and GR 8 exhibit 2D turbulence for the entire U range (Table 16.21) 
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a 




Figure 5.3: The slope a of the H I power spectrum plotted against Ssfr, the SFR 
per unit area. The area has been determined from optical images and H I 
images in the left and right panels respectively. The galaxies with 3D and 
2D turbulence are shown using empty and filled circles respectively. Note 
that the SFR for DDO 210, marked in the figure, is only an upper limit. We 
also include the data for NGC 628, a spiral galaxy. It is marked with an 
arrow. The straight lines show the linear correlation that we find (discussed 
in the text) between the slope of the H I power spectrum and the SFR per 
unit area. 

in the measured H I power spectrum. This imposes the upper limit of 51 pc and 30 pc 
respectively on the scale-height. 

The energy input from st ar formation is belie v ed to be a major driving force for 



the turbulence in the ISM (lElmegreen and Scalo 



20041 ). Also first stage of the star 



formation process is known to be driven by the ISM turbulence. Hence, it is very 
interesting to check whether the slope of the power spectrum of intensity fluctuations in 
these galaxies has any correlation with the SFR. Any correlation, if present, will provide 
an ins ight into the relation between the star formation and ISM turbulence. In a recent 



study 



Willett et al. 



(120051 ) calculates the SFR per unit area for 9 irregular galaxies and 
investigates the correlation with the V-band and Ha power spectra. The length-scales 
they probe are 10 — 400 pc. In Ha they find that the power spectra becomes steeper 



as the SFR per unit area increases. However they do not find any correlation in the 
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V-band. In this paper we probe length-scales 100 pc to 3.5 kpc. We test for a possible 
correlation between the slope of the power spectrum and the SFR per unit area and per 
unit H I mass for the 5 dwarf galaxies in our sample (Figure 15.31) . We also include 
the spiral galaxy NGC 628. We report results using the area estimated from both, the 
H I disk and the optical disk. The data for the SFR, angula r ext ent, inclination in H I 



and H I mass for these galaxies are from 



Begum et al. 



toom and 



Begum et al. 



(2008al). 



The SFR is determined from Ha emission, and the angular extent in H I is determined 
from the column density maps at the value of 10^^ atoms cm ~^. For the angular ex tent 



and inclinatio n of th e optical disk we use the parameters from 



Sharina et al. 



Begum et al. 



f l2006bh and 



(120081 ). The linear correlation coefficient is found to have values 0.35 and 
0.34 for the H I and optical disks respectively, indicating the absence of any correlation 
between the slope of the power spectrum and the SFR per unit area. We have similar 
result for SFR per unit H I mass for these 5 galaxies. 

Turbulence can also possibly be related to different parameters of the galaxy. We 
test for correlations between the slope of the H I power spectrum and the following 
parameters: total H I mass, H I mass to light ratio, total dynamical mass, total baryonic 
mass , gas fraction, baryon fraction. We use the estimates of these observable quantities 



from 



Begum et al. 



( l2008al ). In all the cases the linear correlation coefficient is found to 
lie between —0.5 and 0.5, indicating the absence of correlation. 

In the above analysis we have considered galaxies with both 3D and 2D turbulence 
taken together, which in turn can suppress the correlations. Hence, we further investi- 
gated the correlations by considering the 3 galaxies in our sample with 2D turbulence. 
In this case the correlation coefficient comes out to be —0.70 and —0.74 for SFR per unit 
area of optical and H I di sks respective l y. Th is indicates a strong correlation and the 



result is similar to that of 



Willett et al. 



(120051 ) ■ namely that the surface density of star 



formation rate is larger for galaxies with steeper H I power spectrum. These findings 
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indicate a possible link between star formation and the nature of turbulence in the ISM. 
However, it would not be reahstic to speculate on a possible cause and effect relation 
between these two. It is quite possible that the observed correlation is an outcome of 
extraneous factors which influence both. 



5.4 Summary and conclusions 

In this chapter we have performed and presented the power spectrum analysis of 7 faint 
dwarf galaxies. For 5 of these galaxies we can conclusively estimate the power spectrum. 
Below we summarize the main points of this chapter. 

• ISM of the dwarf galaxies have a scale-invariant power spectrum consistent with 
turbulence. Power law power spectrum of slope —1.1 to —2.6 provides good model 
to these power spectrum. 

• We have observed dichotomy in the slope of the power spectrum which we have 
explained by 3D and 2D turbulence. Based on this we could also provide limits to 
the scale-height of these galaxies. 

• We have investigated possible correlation between the power spectrum slope and 
the dynamical parameters of the galaxies. Our null result suggests that these 
parameters does not influence turbulence in the dwarf galaxies. 

• We have found a correlation between the SFR per unit area (of the optical disk and 
the H I disk) with the power spectrum slope for the galaxies with 2D turbulence in 
our sample. This is a possible indication of the relation between the star formation 
and turbulence in the ISM. However, the Ha emission i s not a good tracer of SFR 



in low mass galaxies due to stochastic star formation f lKarachentsev and Kaisin 
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20071 : iLee et al.l . 120071 ) . Also, the length-scales across which the power spectrum 
has been estimated for the dwarf galaxies in our sample are substantially larger 
than the optical disk where star formation occurs. These points must be taken 
into consideration while interpreting our results. 



Finally, we note that the total number of galaxies in our analysis is rather small for 
a statistical conclusion. A larger galaxy sample is required for a better understanding 
of the generic features, if any, of turbulence in the ISM of faint dwarf galaxies. We 
considered the possibility of performing the power spectrum analysis of a larger sample 
derived from the FIGGS B We tried to estimate the power spectrum of the galaxies 
DDO 43, DDO 125, DDO 189, UGC 7605 and UGC 8508. We found that the FIGGS data 
are noise limited for these galaxies which suggests the necessity of deeper observations. 



* Faint Irregular Galaxies GMRT Survey (iBeeum et al.l . l2008a[) 
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6 Power Spectrum Analysis of the 
THINGS Galaxies 



6.1 Introduction 



In Chapter [3] we have performed power spectrum analysis of three spiral galaxies and 
have found that the power spectrum shows the existence of scale-invariant structures 
in the ISM at a wide range of length-scales starting from ~ 10 kpc to lower. It will 
be interesting to probe ISM of a sample of spiral galaxies to investigate the nature of 
turbulence therein and possible relation of turbulence with the other dynamical proper- 
ties of the galaxy. This requires systematic observations of a hand full of nearby spiral 
galaxies with good spatial resolution. 

THINGS, The H I Nearby Galaxy Surve y, is a survey using H I 21-cm emission of 



20081 ). The main aim of this 



34 nearby galaxies obtained using the VLA (IWalter et al. . 
survey is to investigate the nature of the ISM, galaxy morphology, star formation and 
mass distribution etc. across the Hubble sequence. THINGS provides homogeneous, high 



quality 21-cm data cubes, with high spatial (~ 6) and velocity (< 5.2 km s ^ ) reso 



for a representative sample of local spiral and dwarf galaxies. 



de Blok et al. 



(2008 



ution, 
) have 



75 



6 Power Spectrum Analysis of the THINGS Galaxies 



presented high resolution rotation curve for the galaxies in the THINGS sample. Star 

the T HIGNS galaxies are also very extensively studied 



formation rate and efficiency o 


the T 


(Bigiel et al.. 


2008; 


Lerov et al. 


2008) 



20081 ). This is an ideal data set for a comparative study 



of turbulence in the neutral ISM of spiral galaxies. In this chapter we present the result 
of power spectrum analysis of a sample of 18 spiral galaxies drawn from the THINGS 
survey. 



6.2 Data and analysis 



Here we restrict our analysis to the spiral galaxies in the THINGS sample with minor 
axis greater than 6'. Table 16.11 gives different parameters of the galaxy sample ana- 
lyzed here. The values for the parameters li k e the distance to the galaxy, SFR and the 



total H I mass is taken from 
angle are taken from 



Walter et al. 



Ide Blok et al.l ( 120081 



(12008 



), whereas values for the inclination 



de Blok et al 



(|2008[ ) also present rotation 



curves for the 13 galaxies in our sample. We use these rotation curves to estimate the 



dynamical mass fo r those galaxies. For the rest o: 



ues are taken from 



Bottinelli and Gouguenheim 



Huchtmeier and SeiradakisI (1198 



3 



the g alaxies the dynamical mass va l 



I97I: 



luchtmeier anc 



Kamphuis and Briggd (1199 



3; 



Sofud (Il996[ ) and then 



Witzell (Il979h 



rescaled for the adopted distances noted in Table 16.11 We have estimated the H I 
velocity dispersion using a method described in Section 16.31 Note that, although the 
galaxy NGC 4826 fits our selection criteria, we have not used it in our analysis. It has 
a very bright H I core which makes the window function complicated and the power 
spectrum estimation is not straight forward. 

The observations and data reduction of the THINGS galaxies are discussed in[ 



Walter et al, 



( 12008 



For our analysis, we start with the calibrated uv data prior to continuum sub- 



*We are indebted to Fabian Walter for providing us with calibrated H i data from the THINGS survey. 
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Galaxy 



Major Minor D 



SFR 



HI 



M, 



dy 



ay 





( 


) 


( 


) 


(Mpc) 


( 


) 


(M0 yr- 




(l08 Mq) 


(lO" Mq) 


[km s~ 




NGC 628 


22 





20 





7.3 


15 





1 


2 


38.8 


6.3 


7 


2 


NGC 925 


16 





10 





9.2 


50 





1 


1 


45.8 


1.7 


9 


4 


NGC 2403 


25 





22 





3.2 


55 








9 


25.8 


3.1 


9 


6 


NGC 2841 


22 





26 





14.1 


69 








2 


85.8 


31.2 


21 


6 


NGC 2903 


25 





15 





8.9 


66 









43.5 


4.3 


15 


8 


NGC 3031 


38 





24 





3.6 


59 





1 


1 


36.4 


5.2 


12 


9 


NGC 3184 


17 





14 





11.1 


29 





1 


4 


30.7 


6.4 


9 


1 


NGC 3198 


22 





7 


5 


13.8 
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9 


101.7 


8.1 


11 


6 


NGC 3521 


22 





7 


5 


10.7 


69 





8 


4 


80.2 


12.2 


29 


8 


NGC 3621 


25 





15 





6.6 


62 





2 


1 


70.7 


12.8 


22 


6 


NGC 4736 


18 





12 





4.7 


44 








4 


4.0 


1.2 


11 


4 


NGC 5055 


30 





25 





10.1 


51 





2 


4 


91.0 


12.8 


12 





NGC 5194 


16 





12 





8.0 


30 





6 


1 


25.4 


2.7 


14 


4 


NGC 5236 


30 





24 





4.5 


31 





2 


5 


17.0 


2.6 


10 


9 


NGC 5457 


30 





25 





7.4 


30 





2 


5 


141.7 


5.9 


15 


6 


NGC 6946 


35 





25 





5.9 


35 





4 


8 


41.6 


7.3 


9 


3 


NGC 7793 


12 





9 





3.9 


43 








5 


8.9 


0.7 


9 


9 


IC 2574 


14 





8 





4.0 


51 








1 


14.8 


0.5 


9 


9 



Table 6.1: Some parameters of the galaxies used for the power spectrum analysis. 

Columns 1-9 gives 1) Name of the galaxy, 2) and 3) Major and Minor 
axis at a column density of 10^^ atoms cm~^, 4) Distance to the galaxy, 
5) average HI inclination angle, 6) Star Formation Rate, 7) HI mass, 8) 
Dynamical mass, 9) average velocity dispersion, evaluated from the Mo- 
ment II map. These values are obtained from the following refferences: 



Walter et al. 


f 


20081). 


de Blok et al. 


(2008). 


Bottinelli and Goueuenheim 


(1973r). Huchtmeier and Witzel 


( 


1979) 




Huchtmeier and Seiradakis (19851). 


Kamphuis and Briegs ( 


1992). 


Sofuel ( 


1996). 
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traction. In the standard THINGS pipeline, the continuum is subtracted by fitting a 
linear polynomial to each visibility (i.e. AIPS task UVLIN), whereas for our analysis it 
is important that the strong continuum sources be modelled and then subtracted from 
the visibilities. For each galaxy, the frequency channels with H I emission are identified 
and a continuum image is made by combining all the line-free channels. This contin- 
uum is then subtracted from the data in the uv plane using the AIPS task UVSUB. 
The resulting continuum-subtracted data are used subsequently. We follow the pipeline 
discussed in Section 12.31 for the rest of the analysis. 

6.3 Result and discussion 

For each galaxy, we have identified a baseline range {Umin to Umax) over which the real 
part of V{U) estimated from the channels with H I emission is large compared to both 
its imaginary part and the real part of V{U) estimated from the line-free channels. A 
power law of the form A U°' is found to provide a good fit to the observed power spectra. 
Figure 16. 1[ 16.21 and 16.31 shows the result of our power spectrum analysis. The results 
are also summarized in Table 16. 2^ where Column (2) to (7) gives the value for the 
velocity channel thickness, range of baseline for the fit, corresponding range of length- 
scales and the power law index a respectively for = n (see Section 12.31) . Column (8) 
gives the power law index a for A^ = 1. Unless otherwise stated we will use the values of 
a with N = n for further analysis. The presence of a scale-invariant, power law power 
spectrum indicates that turbulence is operational in the ISM of these galaxies. 

For the galaxies NGC 2841, NGC 3031 and NGC 3521, the angular extent, ^o, is 
comparable to the telescope field of view and hence two observations with two different 
pointing centers (North and South) were used. We refer these two different directions 
using N or S after the galaxy's name. For NGC 3031S, the power spectrum could be 
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Galaxy 




A V 

{km s~^) 


Umin 

(k A) 


Umax 

(k A) 


Rmin 

(kpc) 


Rmax 

(kpc) 


a 

(n = N) 


a 

(n = l) 


NGC 628 




41.6 


1.0 


10.0 


0.8 


7.5 


-1.6 ±0.1 


-1.6 ±0.1 


NGC 925 




41.6 


1.0 


10.0 


0.9 


9.2 


-1.0 ±0.2 


-1.0 ±0.2 


NGC 2403 


* 


83.2 


0.7 


7.0 


0.6 


4.0 


-1.1 ±0.1 


-2.1 ±0.3 


NGC 2841N 




83.2 


1.0 


10.0 


1.4 


14.0 


-1.7±0.2 


-1.7±0.1 


NGC 2841S 




83.2 


1.0 


10.0 


1.4 


14.0 


-1.5 ±0.2 


-1.5 ±0.1 


NGC 2903 




83.2 


0.8 


8.0 


1.1 


11.1 


-1.5 ±0.2 


-1.5 ±0.1 


NGC 3031N 




41.6 


2.0 


10.0 


0.4 


1.8 


-0.7±0.1 


-0.7 ±0.1 


NGC 3184 




41.6 


0.7 


7.0 


1.6 


15.8 


-1.3 ±0.2 


-1.3 ±0.1 


NGC 3198 


* 


83.2 


1.6 


10.0 


1.4 


8.6 


-0.4 ±0.3 


-1.6 ±0.3 


NGC 3521N 




83.2 


1.0 


17.0 


0.6 


10.7 


-1.8 ±0.1 


-1.8 ±0.1 


NGC 3521S 




83.2 


1.0 


17.0 


0.6 


10.7 


-1.6 ±0.2 


-1.6 ±0.2 


NGC 3621 


* 


83.2 


1.0 


12.0 


0.6 


6.6 


-0.8 ±0.2 


-1.9 ±0.2 


NGC 4736 


* 


83.2 


0.6 


10.0 


0.5 


7.8 


-0.3 ±0.2 


-0.9 ±0.2 


NGC 5055 




83.2 


1.0 


10.0 


1.0 


10.0 


-1.6 ±0.1 


-1.6 ±0.1 


NGC 5194 




83.2 


1.0 


8.0 


1.0 


8.0 


-1.7±0.2 


-1.7 ±0.2 


NGC 5236 




83.2 


0.6 


6.0 


0.8 


7.5 


-1.9 ±0.2 


-1.9 ±0.2 


NGC 5457 




83.2 


0.6 


12.0 


0.6 


12.3 


-2.2 ±0.1 


-2.2 ±0.1 


NGC 6946 




20.8 


1.5 


10.0 


0.3 


4.0 


-1.6 ±0.1 


-1.6 ±0.1 


NGC 7793 




41.6 


0.6 


6.0 


0.6 


6.5 


-1.7±0.2 


-1.7±0.1 


IC 2574 




41.6 


1.8 


10.0 


0.4 


3.3 


-1.7±0.3 


-2.0 ±0.2 



Table 6.2: Result of the power spectrum analysis with N = n and = 1. Column 1 to 
8 gives 1) name of the galaxy, 2) width of the channel used to estimate the 
power spectrum, 3) and 4) the range of U value for which the power spectrum 
is evaluated, 5) and 6) correspoindin length scales and 7) the power law index 
a with 1 — 0" error for N = n and Column 8) gives the power law index a 
with 1 — 0" error for N = n. We found the power law index a changes with 
channel thickness for the galaxies marked by (*) in the table. The values for 
the Umax, Umin, Rmax and Rmin are given for = n only 
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estimated over a very limited baseline range and hence we do not consider it here. 

The length-scales over which we could fit the power law (Table 16. 2p covers in total 
a range from 300 pc to 15.8 kpc across the entire sample. The power spectrum of the 
galaxy NGC 5457 is found to be a power law for the length-scale range of 0.6 kpc to 
12.0 kpc which is the largest dynamical range obtained for an individual galaxy in our 
analysis. Note that, NGC 3184 exhibits a power law power spectrum up to a length- 
scale of ~ 16 kpc, where the same for the galaxy NGC 2841 is 14 kpc. These could be 
understood as the existence of turbulence in the ISM of the spiral galaxies at a scale 
comparable to the radius of the galaxy's disk. 

Figure 16.41 shows a histogram of the distribution of the estimated power law index 
a. The histogram is sharply peaked in the range a = —1.8 to —1.5. The mean and 
the standard deviations are —1.4 and 0.5 respectively. We do not find a break in the 
power spectrum of any of the galaxies in our sample. For all the galaxies in our sample, 
the largest length-scales over which the power law is estimated are larger than the 
typical scale-height of galaxies. Following the results of our simulation discussed in the 
Chapter 31 we conclude that we are dealing with 2D turbulence for all these galaxies. 

The galaxies in our sample have a wide range of inclination angles 'i', ranging from 7° 
for NGC 628 to 74° for NGC 2841. To investigate whether the spread in the value of a 
seen here is an effect of the difference in i, we evaluate the linear correlation between the 



power law index a and the inclination angle of the galaxy (Figure 16. Sp . The resultant 
linear correlation coefficient of 0.25 is a clear indication that the dispersion in the values 
of a is not due to the inclination angles. 

At small scales, a major driving force of the ISM turbulenc e is the energ y inpu t from 



star formation processes. For a sample of irregular galaxies 



Willett et al. 



(120051) finds 



that on scales of 10 — 400 pc the power law index of the Ha power spectra becomes 
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Figure 6.1: The real (triangle) and imaginary (circle) part of the estimator evaluated 
using the channels with HI emission is plotted with the 1 a error bars in a 
log log scale against the baseline values for the galaxies NGC 628, NGC 925, 
NGC 2403, NGC 2903, NGC 2841, NGC 3031 and NGC 3184. The real 
part evaluated with the line-free channels (star) is also shown. The best 
fit power law is given by the bold solid line. Note that in the U range the 
power spectrum is deduced, the real part of the estimator has a higher value 
compared to the imaginary part or the real part of it from the line free 
channels. 
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Figure 6.2: The real (triangle) and imaginary (circle) part of the estimator evaluated 
using the channels with HI emission is plotted with the 1 a error bars in log- 
log scale against the baseline values for the galaxies NGC 3521, NGC 3198, 
NGC 3621, NGC 4736, NGC 5055, NGC 5194 and NGC 5236. The real 
part evaluated with the line-free channels (star) is also shown. The best 
fit power law is given by the bold solid line. Note that in the U range the 
power spectrum is deduced, the real part of the estimator has a higher value 
compared to the imaginary part as well as the real part of it from the line 
free channels. 
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U [kAl 



Figure 6.3: The real (triangle) and imaginary (circle) part of the estimator evaluated 
using the channels with HI emission is plotted with the 1 a error bars in log- 
log scale against the baseline values for the galaxies NGC 5457, NGC 6946, 
NGC 7793 and IC 2574. The real part evaluated with the line free channels 
(star) is also shown. The best fit power law is given by the bold solid line. 
Note that in the U range the power spectrum is deduced, the real part of 
the estimator has a higher value compared to the imaginary part as well as 
the real part of it from the line free channels. 
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Figure 6.4: Shaded regions shows the histogram of the distribution of the power law 
index a. Note that the histogram has 7 bins for 18 data points. Sample 
mean fi and sample standard deviation a is also shown. 



r = 0.25 



20 



40 
( ° ) 



60 



Figure 6.5: Scatter plot of the average inclination angle i with power law index a. The 
1 a error bars of a are also shown. The value of the linear correlation 
coefficient r is given at the top right corner. 
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r = -0.07 



-6 
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Figure 6.6: Scatter plot of the surface density of the star formation rate (Table 16. ip 
with power law index a. The 1 a error bars of a are also shown. The left 
and right pannels shows the SFR per unit area of the optical and HI disk of 
the galaxy respectively. The value of the Linear Correlation Coefficient r is 
given at the top right corner of the each panel. 

steeper as the SFR per unit area increases. In Chapter O we have reported a weak 
correlation between the SFR and the power law index of the H I intensity fluctuations 
in dwarf galaxies. We have calculated the linear correlation coefficient of a with the SFR 
per unit area of the optical disk as well as the SFR per unit area of the H I disk. The 



result, shown in Figure [6. 6[ suggests that on large-scales at least, there is no correlation 
between the SFR and the power law spectral index. 

We next investigate the linear correlation of the power law index a with the large scale 
dynamical parameters of the galaxies namely the H I mass {Mm) and the dynamical 
mass (Mdy). We use the values of the maximum rotatioi ial velocity (Vmnx) a-i id the maxi- 



mum radius (rmax) from the rotation curves presented in 
the dynamica l mass enclosed within that radius, viz. {M^y 



de Blok et al. 



de Blok et al. 



*^ max ' max 



(120081) to evaluate 
/G). Note that 



(j2008h have given rotation curves for only 13 galaxies in our sample. For 



the other 5 galaxies, namely NGC 68, NGC 3184, NGC 5194 NGC 5236 and NGC 5457 



we have taken the M^„ estimates from earlier references, viz. (IBottinel 



1973 



Huchtmeier and Witzel 



1979 



Huchtmeier and Seiradakis 



1985 



i and Gouguenheiml. 



Kamphuis and Briggs 
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Figure 6.7: Scatter plot of the dynamical mass M^y and the total HI mass M^i with 
power law index a. The 1 a error bars of a are also shown. The value of the 
Linear correlation coefficient r is given at the top right corner . 

' The Mhi values are taken from 



1992 



Sofue 



199( 



Walter et al. 



mm . Figure Em 



shows the scatter plot of a with M^y (left) and Mhi (right) respectively. Our analysis 
shows that a is not correlated with the dynamical mass or the total H I mass. 

ISM turbulence can also be driven by the kinetic energy of the H I gas. In this case, 
one would expec t that the H I veloci ty dispersion (a^j) would be correlated with a. We 



use Eqn. (7) in 



Walter et al 



(2008 



) to define the average H I velocity dispersion. In 



estimating ay, we perform a 8 pixel smoothing on the THINGS archival Robust Weight 
velocity dispersion (MOMENT 2) maps HI and then resample it to a resolution of 6". The 
overall velocity dispersion of the galaxies in the sample is evaluated using the values in 
individual pixels of the resampled map (Table 16. 1^ Column (9). Figure 16.81 shows the 
scatter plot of the cr„ values with a. We find no correlation between a and a^. 

In summary, the slope a is not correlated with any of the galaxy parameters that we 
have considered. While ~ 40% of the a values are in the narrow range —1.8 to —1.5, 
the remaining galaxies have a values spread over the broad range —2.2 to —0.3. It is 
unclear if this arises from statistical fiuctuations in the distribution of a, or if there is a 
physical mechanism underlying this spread. 



■|'http://www.mpia-hd.mpg.dc/THINGS/Data.html 
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Figure 6.8: Scatter plot of the average velocity dispersion with power law index a. 

The 1 a error bars of a are also shown. The value of the Linear correlation 
coefficient r is given at the top right corner. 



Both H I density fluctuations as well as spatial fluctuations in the ve 
gas co uld contribute to observed fluctuations in the H I specific intensity. 



ocity of the H I 



Lazarian and Pogosvan 



(120001 ) have shown that it is possible to disentangle these two contributions by studying 
the behavior of the H I power spectrum as the thickness of the frequency channel is 
varied. If the channel thickness Av is greater than the velocity dispersion for a galaxy, 
the velocity fiuctuations are averaged out and the intensity fiuctuation power spectrum 
becomes a direct probe of the density fiuctuation alone. We investigate the possible 
effect of the velocity fiuctuation in the observed intensity fiuctuation power spectrum by 
estimating Phi{U) increasing the channel thickness (Av) beyond the velocity dispersion 
of the galaxy. We find the slope a of the power spectrum of the galaxies NGC 2403, 
NGC 3198, NGC 3621 and NGC 4736 increase as the channel thickness is varied beyond 
a value Av = v^. These galaxies are marked by an Astrix in the Table 16. 2[ We show 
the variation of a with the increasing channel thickness for these galaxies in Figure 16. 9[ 
It is interesting to note that for all four galaxies Vc ~ 20km s~^ and the change in slope 
is independent of the velocity dispersion of the galaxy in question. We interpret the 



value of a for the largest Av to be the power law index of the density fiuctuation power 



87 



6 Power Spectrum Analysis of the THINGS Galaxies 
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Figure 6.9: Change in the slope of the power law index a with velocity channel thick- 
ness is shown for the four galaxies NGC 2403, NGC 3198, NGC 3621 and 
NGC 4736. Levels at the bottom shows the width of the velocity channels 
whereas the top levels show the number of channels averaged together before 
evaluating the power spectrum. 
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spectrum (Figure 16.91) . The slope of the power law is found to remain constant as the 
thickness of the channels is varied for the rest of the galaxies in the sample. Since for 
all the galaxies the thickest channel is considerably larger than the velocity dispersion, 
we conclude that for these galaxies, the H I power spectrum is purely due to density 
fluctuations and not the velocity fluctuations in the H I gas. 

6.4 Summary and conclusion 

This paper summarizes the power spectrum analysis of 18 spiral galaxies from the 
THINGS sample. The estimated power spectra can be well fitted with a power law, 
indicating the presence of turbulence in the ISM. This analysis is the first comprehen- 
sive study of turbulence in spiral galaxies with a moderate sample size. In this section 
we give a brief summary of the important findings from our power spectrum analysis. 

• We find that the galaxies in our sample shows a power law power spectrum over 
length-scales ranging from ~ 300 pc to ~ 16 k pc with several galaxies having 
Rmax > 10 kpc. Presently there is no known mechanism which can drive turbulence 
at such large length-scales. 

• We have observed that for most of the galaxies in our sample, a ~ —1.6 (Fig- 
ure 16. 4p . This value of a is similar to that is found for 2D power spectrum of the 
galaxies in Chapter [31 

• In Chapter O we have found no significant correlation between a and the SFR, 
dynamical mass, total H I mass of dwarf galaxies. Here we estimate the power 
spectrum of 18 external spiral galaxies which provide us a moderate sample size. 
However, our analysis also shows that the a has no correlation with the inclination 
angle, SFR, dynamical and total H I mass of these galaxies. 
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6 Power Spectrum Analysis of the THINGS Galaxies 



For most of the galaxies in our sample the width of the frequency channel is found 
to have no effect on the power law index. The slope of the power spectrum re- 
mains the same as the channel-thickness becomes greater than the average velocity 
dispersion of the galaxies. This suggests that the intensity fluctuation power spec- 
trum is a direct probe of the column- density fluctuations in the galaxy's disk. 
Interestingly, we notice that for four galaxies NGC 2403, NGC 3198, NGC 3621 
and NGC 4736 the slope of the power spectrum decreases as the channel width 
Av is increased, and beyond a certain value of Av ~ 20 km s~^ it remains con- 
stant. However, the channel width where this transition occurs has no correlation 
with the velocity dispersion of th e galax y considered. We found this to be differ- 



ent from 



Lazarian and PogosvanI (120001 ) and presently do not have any physical 



understanding of the reason behind this. 
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7 Summary and Future Scopes 



This thesis presents our effort of probing turbulence in the ISM of the external galaxies 
by estimating the power spectrum using radio-interferometric observations of H I 21-cm 
radiation. Here we summarize our main results and discuss the future scopes. 

• We have estimated the power spectrum of 25 external galaxies including 5 nearby 
dwarf and 20 spirals. For every galaxy, the power spectrum is found to follow a 
power law over a wide range of length-scales which spans 60 pc to 16 kpc in total. 
This indicates that the ISM is turbulent at length-scale comparable to the radius 
of the galaxy's disk. 

• We have observed a dichotomy in the estimated value of the power spectral index 
and have explained it in terms of 2D and 3D turbulence at large and small scales 
in the galaxy's disk respectively. Based on a break in the power spectrum of the 
galaxy NGC 1058 we have estimated its scale-height as 450 ± 60 pc. This is the 
first determination of the scale-height of any external spiral galaxy. 

• We investigated whether the turbulence is related to the different dynamical pa- 
rameters like dynamical and H I mass, star formation rate etc. of the galaxies. 
For the dwarf galaxies we found a weak correlation between surface density of star 
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7 Summary and Future Scopes 

formation rate and the power law index. There was no such correlation observed 
for the spiral galaxies. 

• The observed power spectrum is found to have no effect of the velocity fluctuations 
for 21 of the 25 galaxies we have analyzed. For the rest of the galaxies we have 
observed power law index changes with the channel thickness indicating an effect 
of the velocity fluctuation. 

• Our power spectrum analysis of the harassed galaxy NGC 4254 has shown that the 
galaxy harassment effects the ISM structures. The central part of such a galaxy 
exhibit different dynamics compared to that in the outer disk. 

• Given a power law power spectrum, its estimator introduced in this thesis is sensi- 
tive to the slope but not to the amplitude of the power law. On the other hand, for 
smaller length scales in our galaxy, the amplitude and slope of the power spectrum 
is very well constrained. If one can carry out observations with more sensitivity 
with better uv-coverage at the larger baseline range, it would be possible to probe 
similar length scales in the external nearby galaxies with the existing radio interfer- 
ometers like VLA or GMRT. This will provide a way of calibrating the amplitude 
of the estimator we have used and then a direct comparison between the large and 
small scale powers can be done. Also note that we have been restricted to the 
nearby galaxies because of the limited angular resolution and hence linear resolu- 
tion in the disk of the galaxies that we can achieve with the present telescopes. 
For an even better uv-coverage and more dynamic range in the length scale within 
which the power spectrum can be estimated, upcoming radio telescopes like SKA 
will come be very useful. 
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The known energy sources of ISM turbulence include galaxy rotation, self gravity, 
supernovae etc., which are believed to drive turbulence at parsec scales. It is not clear if 
these are equally effective in explaining the large scale turbulence observations reported 
in this thesis. A detailed physical modeling is required to have better understanding 
of ISM dynamics. It may be possible to estimate the ISM velocity fluctuation power 
spectrum, which can independently test the existence and nature of turbulence in the 
ISM. We wish to perform numerical simulation to understand the effect of velocity 
fluctuations etc. in our power spectrum estimation. Finally, it will also be interesting 
to study the ISM turbulence in a sample of galaxies residing in the cluster environment. 
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A Incomplete uv coverage and correlated noise in the 
radio-interferometric observations 

Noise in the visibility plane J^iU) corresponds to I^^\0) given by 

lW{0) = J dU e'^^^-^MiU) (Al) 

Fourier transform is a linear operation and in practice a discrete Fourier transform is 
used. We can write Eqn. (lAlj) as 

k 

where Fj^ are the components of the transformation matrix Fj^, = exp(^^^). In a 
typical radio-interferometric observation the uv plane is not completely sampled. The 
actual inverse-transformation from the measured visibility to the specific intensity is 
given by 

/f) = Y.kVkiNi, (A3) 

where T^iNi gives the measured visibility points. We call this matrix F as the projection 
matrix. It is to note here, that F is an Identity matrix if all the points in the uv plane 
are sampled, otherwise it is a diagonal matrix with diagonal elements as or 1. In the 
latter case the inverse of this matrix does not exist. Components of the noise covariance 
matrix (assumed to be diagonal) can be written as 

N,, = (A4) 

We calculate the auto-correlation function as 

= VitTtr^lTl^^ij. (A5) 



7 Appendix 



In case of complete uv coverage, we have 



= ^•5.,, (A6) 

since F is unitary. However for real observations with limited uv coverage T^p is not a 
identity matrix. The noise power spectrum is given by, 

mj mi SiJ 

— ^mp F^^- (A7) 

This clearly shows that for limited baseline coverage, there exists correlated noise in the 
image plane. Every attempt to estimate power spectrum from the Fourier transform 
estimators hence has to take care of this correlated noise effect. 
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B Validity of the 2D approximation. 

In Chapter [2] we have assumed the basehne ?7 as a 2D vector. This is rather an 
approximation and in reahty 
vahdity of the 2D assumption, 



approximation and in reahty the vector U has three components. Here we discuss the 



Bl 2D and 3D visibility, the effect of u;-term on the aperture 

The direct observable in the radio-interferometric observations is the complex visibility 
V^^{U, u). For a pair of antennae separated by d, with each antenna pointing along the 
direction of the unit vector k (referred to as the phase center) we have 

V3^(t/,z/) = e^™ j d^lne^^''^-^'^'^^A{n-k,u)I{n-k,p) (Bl) 

where n denotes the unit vector to different directions of the sky, baseline U = d/X, 
A{n — k, v) denotes the primary beam and J(n — fc, v) is the specific intensity. Writing 
U = U± + wk, where U± is a 2D vector, and defining n — k = 6, we have, for | |^ 1, 
6 ■ k ^ 0, implying that is a 2D vector. In this limit 6 gives the position of any 
point on the sky with respect to the phase centre in a 2D tangent plane. This is known 
as the flat-sky approximation. The term wk quantifies deviation from this. In the 2D 
approximation we have 

V^^{U±, u) = j d0 e^^'^^-^ A{e, v) I{e, v) (B2) 

Writing the specific intensity /(0, v) as /(0, z/) = Jj, + 5/(0, v\ where the first term is a 
constant background and the second term is a fluctuation, we have 

V'''(C/±, v) = I,A{U^, u) + i(t/^, u) ® SI{U^, u) (B3) 

where tilde represents a Fourier transform and ® denotes a convolution. 

The aperture function A{U u), peaks at U± = and has a finite width. Hence, we 
shall retain the second term in all subsequent discussions. 



*Thc Appendix is adopted from the originally published in t he paper titled "The Effect of w-term on 
Visibility Correlation and Power Spectrum Estimation" bv lDutta et al. ( 2010b ). Note that we have 



used different symbols here compared to what is defined in rest of the thesis for reasons given in the 
text. 
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Figure Bl: Schematic representation of a particular visibility measurement with a pair 
of antenna A and B separated by a distance \ d \. The square box at the 
top right corner represents the portion of sky observed or the field of view 
(FOV). The black dot at the centre of the FOV corresponds to the direction 
to the phase centre k. Intensity fluctuations at a given direction n to the 
sky is denoted as I{n — k, u) or 1(0, u). Note that is a 2D vector in the 
plane of the sky. 
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We note that, ignoring the w-term leads to a simphfication of the expression for visibil- 
ity and in the 2D approximation, V^^(C/_l, u) is the Fourier transform of A{0, h')5I{0, v). 
Hence, we have, 

A(6», u) 51(6, u) = j dU'^ V2^(t/^, z/) e-^^'^'^-^ (B4) 
Substituting in Eqn. ( IBlj) we obtain 

V3«(C/, u) = j dU'^ K{U, U'^, v) V'^(C/'^, v) (B5) 
Where the kernel K{U , U'j_, v) is given as 

= 47rjo(27r|C/'^-C/|), (B6) 

with Jo denoting the Spherical Bessel function. 
Defining a quantity A(U, u) as 

A{U, z/) = 47r j dU'^ Jo {2n \ U'^ - U |) A{U'^, v). (B7) 

The scalar 3D visibility takes the form 

V^^(C/, v) = j dU'^ A{U - U'^, v) 5I{U'^) (B8) 

It is to be noted that the 'w' dependence of V^^{U,v) is translated to the function 
A{U — t/_L,i^). This can be regarded as a modified aperture function. Note that 

valued function, unlike the complex Gaussian as discussed in 



this function is a rea^ 



Hobson and Maisingerl (120021 ). We investigate the nature of the modified aperture 



A{U ^ v) as a function of U±_ for different values of w. We have assumed the primary 



aperture AiXJ u) to be a Gaussian, exp 



of width Uq, and evaluated the integral 



in Eqn. ( IBTj) numerically. We shall not explicitly write the u dependence of A and A 
henceforth. It can be shown that the modified aperture preserves the circular symmetry 
of the primary aperture. The implicit u dependence is present in these functions through 
Uq. For an antenna of diameter D, Uq can be approximately written as f/o ~ -D/A for 
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V = 150.0 MHz 



Figure B2: Modified aperture A{U ^^vS) plotted as a function of U ^ for two different 
values of w, (A) w = 400. A and (B) w; = 12. kA at z/ = 150 MHz. Solid line 
in B shows the Gaussian envelope. 

observing wavelength A. We use Uq = D /\ with D = 45 m (specifications of the GMRT) 
for the subsequent discussion. This corresponds to a FOV of 4° at 150 MHz. We shall 
discuss the effect of larger FOV later. 

Figure |B2] shows the numerically evaluated modified aperture A{U w) as a function 
of U± for two representative values of w (A: w = 400 X, B: w = 12 kA) at frequency 
1/ = 150 MHz. 

For w = one reproduces the primary Gaussian aperture trivially. We also note that 
for w <^ U± the Gaussian profile is still maintained (Figure IB2[ A). However for large 
values of w the aperture function manifests oscillations (Figure IB2[ B). The period of 
these oscillations is found to be sensitive to w (decreasing as w increases). The envelope 
of the modified aperture is maximum at U± = and falls of gradually as U± increases. 



We found that a Gaussian, C exp 



provides a good fit to the envelope. Ignoring 



the effect of oscillations in A, we note that the 3D formalism can be recast in the same 



form as its 2D counterpart with [/^ taking the role of Uq. iHobson and Maisingerl (j2002[ ) 
have obtained a similar result assuming small FOV, where, they have shown that the 
effect of ty-distortion can be considered as turning the primary beam into a complex 
Gaussian. 

We next investigate the dependence of f/^ on w and frequency. For a given frequency 
at large w, is found to increase linearly with w, i.e, Uyj ~ m{v)w. This implies 
that the ty-term effectively broadens the aperture of the instrument. Figure IB3I shows 
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u 1 , 




u - 


150.0 MHz 






235.0 MHz 




1/ — 


320.0 MHz 




u - 


610.0 MHz 




J/ — 


1420.0 MHz 



^ 1 10 100 

w (k A) 



Figure B3: is plotted as a function of w for different central frequencies of GMRT. 
Horizontal solid line corresponds to the maximum possible baseline. 



o 
o 
o 




200 400 600 

U ^ (A) 



Figure B4: Modified aperture A{U±,w) plotted as a function of U± for w = 120. k A 
at = 1420. MHz. 
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the variation of f/^ with w for different frequencies in log- log scale. We note that, the 
slope m(i/) (as represented by the y-intercept in the Figure IB3|) determines the effect 
of w-teim for increasing values of w. A small value of m(i^) implies a slow increase of 
Uyj with w and the effect of the w-term is less. m(z/) is found to fall off as ~ 1/z^ with 
frequency. This indicates that the departure from the flat-sky approximation is more 
pronounced at the lower frequencies. Redshifted 21-cm line observed at frequency u 
probes the redshift z = ^^^'^(^^^) _ i _ Hence, one may expect the w-term to have a 
greater effect while probing higher red-shifts. 

In this thesis we use 21-cm line to study the ISM dynamics of the nearby galaxies 
(z ~ 0). Figure IB4I shows the modified aperture function for frequency 1420 MHz and 
w = 120 kA (this being the maximum U for GMRT like arrays). At this frequency, 
for GMRT, Uq = O.lkA, whereas \w=u„,ax= O.lSkA. It follows that the fiat-sky 
approximation can be safely used if the largest length scale probed, corresponds to a 

Till now, we have investigated the effect of u^-term using D = 45 m, which corresponds 
to an FOV of 4° at 150 MHz. We estimated U^j \w=Umaa: assuming D = 4 m to D = 45 
m aX u = 150 MHz. At D = A m, (which corresponds to the largest proposed FOV, 45°, 
of SKA), \w=u„iax= 4 kA for Umax = 120 kA. For all the cases mentioned above, a 



Gaussian of the form C exp 



provides a good fit to the envelope of the modified 



aperture function and we could set a value for Uw Further, we also observe that for 
these cases |u;=(7^^^~ ^^^5 where dmax is the largest baseline of the telescope array. 
This implies that Uu, \u]=Umax depends only on the geometry of the array configuration. 



B2 Visibility correlation and power spectrum estimation 

For 3D visibilities, power spectrum P{U±) can be estimated by performing the visibility 
correlation as 

Vi^'iU) = (V=^^(C/)V3^(C/)*) = J d^U'^ \A{U - U'^)f P{U[) (B9) 

Noting that the effect of the w-term is contained in the modified aperture A. We can re- 
trieve the 2D estimator defined in Section 12.21 (consider same baseline correlation) 
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s 

> 



0.01 0.1 1 10 100 

Ui(A) 

Figure B5: as a function of U± for w = (dot-dash) and w = Umax (dash) at 
= 150 MHz, assuming P{U^) = U^'^. We also plot P{U^) = U^'^ (solid 
line) for reference. The vertical arrows show the U± value above which the 
power law is recovered at 1%. Note that the plots are given arbitrary offset 
for clarity. 

by replacing A with A. Hence, 

Vi^'iU^) = I d'u'^ \A{U^ - U'^)f P{U'^). (BIO) 

We note that V^^{U) and ^^^(t/i) become V^^iU) and Vi^{UA_) for circularly sym- 
metric primary aperture. We shall now discuss the effect of the w-term on the estimator 
defined in Eqn. (IB9P . 

Figure IB5I shows plotted as a function of for two values of w, (i/; = and 
w = Umax) at z/ = 150 MHz, assuming P{U±) = f/j_~^. We have chosen Umax = l2kX 
(this being the largest baseline for the GMRT at 150 MHz). We have shown the power 
spectrum P{U±) = U±~'^ for comparison. For large values of w, show oscillations 
for U± < Uwi which arises due to the oscillatory nature of |>l(t/)p. 

We find that Vj^^ faithfully recovers the power law f/_L~^ (at 1%) for U±_ greater than 
a certain value. This value is found to be 3 x (v^f/o) for w = and 3 x [■\/2Uyj) at 
w = 12 kA. Hence, a non-zero w-term changes the Ux_ value beyond which the power 
spectrum estimation would be valid. 

The quantity of interest in power spectrum estimation using the radio-interferometric 
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w = 


W ~ Umax 


Aperture 






width 


Uo 




Visibility 






correlation 







Table Bl: The effect of w;-term, comparison between various quantities. 





150 (M Hz) 


1420 (M Hz) 




0.01 


0.1 


Umax 


12.0 


120.0 


Uw \w = Umax 


0.18 


0.18 



Table B2: Relevant U±_ (k A) values at different frequencies, 
observations used in this thesis is actually the quantity 

rUmax 

£[Ui_)= I dw p{w) Vi^{U±,w), (Bll) 
Jo 

where, p{w) is a normalized probability distribution of w. The function p{w) is specific 
to an observation as well as to the array configuration of the interferometer. Hence, 
it is difficult to make a general quantitative statement regarding the effect of w-term 
in S. Since, for a given w, the largest baseline above which V2^{U±,w) ~ P{U±) is 
Uw, we can qualitatively state that S gives a good estimation of the power spectrum 
for ?7 > 3 Uw \w=Umax- It is important to note that, for a specific array configuration, 
Uw \w=Umax is independent of the frequency u, whereas Umax v (Table IBip . Hence, 
the U ^ range amenable for power spectrum estimation is larger for large observing 
frequencies and the use of the 2D visibility for the analysis presented in this thesis is 
justified. 
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C Error in the visibility correlation estimator 

The real and imaginary parts of the estimator V{U) have uncertainties arising from (1.) 
the sample variance and (2.) the system noise. In this section we provide a method of 
estimating the error in the visibility correlation estimator. 



CI The noise covariance matrix 

In radio astronomy the total power available at a radio telescope terminal is conveniently 
stated in terms of the system temperature T'*^* = (Total power )//cb, where is the 
Boltzmann's constant. The system temperature obtained for the antenna looking at a 
blank sky gives a measure of the total random noise in the system. This noise is assumed 
to be Gaussian random with zero mean and its variance for the j^^ antenna can be given 

by a] = T;y\ 

It can be shown that for most practical purposes, a pair of baselines i and j observing 
a point source with flux densit y S, with an integration time of T, at a bandwidth of Au 



the noise covariance matrix is ( IChengalur et al.l . 120071 ) 



rpsys rpsys 

It is important to remember that usually the interferometers do not record the self 
correlations and hence all i = j are not considered from the above expressions. Now 
if the antenna pairs i — j corresponds to a baseline of U, then variance of both the 
real and imaginary part of N'{U,v) is 0"^. For a typical observation, the noise is also 
uncorrelated and constant across the frequency channels. 

An important property of the noise is that at two baselines U a and U b the noise is 
uncorrelated, i.e, 

{M{Ua)M{Ub)) = 5ab N2(C/a) (C2) 

where N2(C/yi) is noise covariance matrix for a given baseline Ua- Assuming Poisson 
system noise we can write. 



since N2(C/) is independent of U. 
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C2 Variance in the real and imaginary part of V{U) 



We assume that the density fluctuations produced due to the ISM turbulence is ho- 
mogeneous and isotropic. In such cases we can replace the ensemble average in the 
Eqn. (I2.12j) by an angular average. The real part of each multiplication of the visibili- 
ties V{Ui) and V{Ui + AU) with AU < l/n9o, provide us an independent estimate, 

V{Ui) = Re [ViUi)V*iUi + AU)] = Vi. (C4) 

We consider annular regions in the uv plane as shown in Figure IC2I (the region between 
Ua and Ub) and average over all the individual estimates Vi in the annulus to get 

1=1 

where, Ni, is the number of individual estimates Vi in the given U bin. This value of V 
is then proportional to the power spectrum Phj{U) for U = {Ua + Ub)/2. 

We will now calculate the variance of the power spectrum estimator we have just 
discussed. It is to be noted that, since visibilities at two different baselines Ui and 
Ui + AU remains correlated for AU < I/ttOq, all the estimates Vi of Eqn. (IC5p are 



not independent. Figure ICll shows three different estimates of V, namely Vi, V2 and 
V3, among which Vi, V3 and V2, V3 are mutually independent estimates, but, Vi, V2 
are not independent. This suggests that we can break the Pj's in say Ng groups such 
that, inside a group, the "Pj's are not independent, however, the estimates from different 
groups are independent. Hence, we can write, 

1 

9 ■ 1 

1=1 

where, Ni are the number of P's in the i*^ group and 

^^^ = W-^^^ ^^^^ 

with {Vn^Vn,) = for z 7^ j and J^'^J^ = ^6- 
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= lie %") 




= lie (V^ 







Figure CI: A particular annular region between Ua and Ub is shown. Note that the 
visibility V{Ua) and V{Un) are correlated with all the visibilities in the 
corresponding shaded circular regions. The estimates Vi and V2 are not 
mutually independent estimates where as the estimates Vi and V3 or V2 and 
V3 are mutually independent estimates. 
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Variance of the estimator is defined as 

where a^^ = {V%^) - {Vn^^. Assuming V2(Af/) = {V{Ui)V{Ui + 5U)), we can write, 

= 2 + 2iV^ ^^^^ 

Hence, using Eqn. ( ICSj) . we get the variance of the estimator given by 



4 = {v^)-{v? 

(CIO) 



V2(AC/)2 + V2(0)2 ^ Nl + 2V2{0)N2 



2N, 2iV2 ^ 

In general, all A^j's are different. If we assume Ni = N2 = ■ ■ ■ = Aj, then above 
Eqn. ( IClOj) becomes 



V2 {A Uf + V2{0y 2V2(0)A^2 + Nl 

2Ng ^ N, 



^ .2V^^; J.2Vw; ^..2lw;.v2^iV2 ^^^^^ 



For an observation, when A'2 ^ V2(0), we can approximate the above equation as 

, V2(At/)^ + V2(0)^ , ATI 
= 2A^ + 2N, ^^'^^ 

Since, we expect the value of the correlation V2(AU") at two different baselines separated 
by I 5U \« is approximately same as the correlation V2(0) at the same baselines, 
we can write. 



This expression is used by lAli et al.l (|2008[ ) in calculating the variance in the power 



spectrum of the foreground signal in EOR (Epoch of Reionization) observation. We can 
identify the first term in Eqn. (1C13[) as a contribution from the sample variance and 
the second term arising due to the statistical fluctuation. 

For the sake of completeness we mention here that the visibility correlation defined in 



120 



C Error in the visibility correlation estimator 



Eqn. (I2.13|) also have an imaginary part defined as 

V'iUi) = Im [V{U^)V*{U, + AU)] = Vl. (C14) 
However, when averaged over all the realizations, it can be shown that, 

(^') = (iv E^/) = o (CIS) 

with 

Since, we expect the value of the correlation V2(AC/) at two different baselines separated 
by I 5U |<< ^Q^^ is approximately same as the correlation V2(0) at the same baselines, 
we can safely neglect the first term compared to the second term. Further, if Ni = N2 = 
■■■ = Ni, and for N2 > 52(0), Eqn. (ICT6D becomes 

4, - ^ (C17) 
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